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Abstract —In this paper we study a family of controllers that 
guarantees attitude synchronization for a network of agents 
in the unit sphere domain, i.e., S 2 . We propose distributed 
continuous controllers for elements whose dynamics are con¬ 
trollable, i.e., control with torque as command, and which can 
be implemented by each individual agent without the need 
of a common global orientation frame among the network, 
i.e., it requires only local information that can be measured 
by each individual agent from its own orientation frame. The 
controllers are constructed as functions of distance functions 
in S 2 , and we provide conditions on those distance functions 
that guarantee that i) a synchronized network of agents is 
locally asymptotically stable for an arbitrary connected network 
graph; ii) a synchronized network is asymptotically achieved 
for almost all initial conditions in a tree network graph. When 
performing synchronization along a principal axis, we propose 
controllers that do not require full torque, but rather torque 
orthogonal to that principal axis; while for synchronization 
along other axes, the proposed controllers require full torque. 
We also study the equilibria configurations that come with 
specific types of network graphs. The proposed strategies can be 
used in attitude synchronization of swarms of under actuated 
rigid bodies, such as satellites. 

I. Introduction 

Decentralized control in a multi-agent environment has 
been a topic of active research for the last decade, with appli¬ 
cations in large scale robotic systems. Attitude synchroniza¬ 
tion in satellite formations is one of those applications [1], 
where the control goal is to guarantee that a network of 
fully actuated rigid bodies acquires a common attitude. 
Coordination of underwater vehicles in ocean exploration 
missions can also be casted as an attitude synchronization 
problem [2], 

In the literature of attitude synchronization, different so¬ 
lutions for consensus in the special orthogonal group are 
found [1], [3]—[10], which focus on complete attitude syn¬ 
chronization. In this paper, we focus on incomplete attitude 
synchronization, which has not received the same attention: 
in this scenario, each rigid body has a main direction and the 
global objective is to guarantee alignment of all rigid bodies’ 
main directions; the space orthogonal to each main direction 
can be left free of actuation or controlled to accomplish 
some other goals. Complete attitude synchronization requires 
more measurements when compared to incomplete attitude 
synchronization, and it might be the case that a rigid body 
is not fully actuated but rather only actuated in the space 
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orthogonal to a specific direction, in which case incomplete 
attitude synchronization is still feasible. Incomplete atti¬ 
tude synchronization is also denoted synchronization on the 
sphere in [11]—[16], where the focus has been on kinematic 
or point mass dynamic agents, i.e., dynamical agents without 
moment of inertia. 

In [4], attitude control in a leader-follower network of rigid 
bodies has been studied, with the special orthogonal group 
being parametrized with Modified Rodrigues Parameters. 
The proposed solution guarantees attitude synchronization 
for connected graphs, but it requires all rigid bodies to be 
aware of a common and global orientation frame. In [5], 
[6], a controller for a single-leader single-follower network 
is proposed that guarantees global attitude synchronization 
at the cost of introducing a discontinuity in the control 
laws. In [7], attitude synchronization in a leader-follower net¬ 
work is accomplished by designing a non-linear distributed 
observer for the leader. In [17], [18], a combination of 
a tracking input and a synchronization input is used; the 
tracking input adds robustness if connectivity is lost and 
it is designed in the spirit of leader-following, where the 
leader is a virtual one and it encapsulates a desired trajectory; 
however, this strategy requires all agents to be aware of 
a common and global reference frame. In another line of 
work, in [3], [8], attitude synchronization is accomplished 
without the need of a common orientation frame among 
agents. Additionally, in [3], a controller for switching and 
directed network topologies is proposed, and local stability 
of consensus in connected graphs is guaranteed, provided 
that the control gain is sufficiently high. In [1], attitude 
synchronization is accomplished with controllers based on 
behavior based approaches and for a bidirectional ring 
topology. The special orthogonal group is parametrized with 
quaternions, and the proposed strategy also requires a com¬ 
mon attitude frame among agents. In [19], a quaternion 
based controller is proposed that guarantees a synchronized 
network of rigid bodies is a global equilibrium configuration, 
provided that the network graph is acyclic. This comes at the 
cost of having to design discontinuous (hybrid) controllers. 
A discrete time protocol for complete synchronization of 
kinematic agents is found in [20]. The authors introduce 
the notion of reshaping function , and a similar concept is 
presented in this manuscript. The protocol provides almost 
global convergence to a synchronized configuration, which 
relies on proving that all other equilibria configurations, 
apart from the equilibria configuration where agents are 
synchronized, are unstable. In [9], controllers for complete 
attitude synchronization and for switching topologies are 
proposed, but this is accomplished at the kinematic level. 



i.e., by controlling the agents’ angular velocity (rather than 
their torque). This work is extended in [10] by providing 
controllers at the torque level, and similarly to [1], stability 
properties rely of high gain controllers. 

In [11], [12], incomplete synchronization of kinematic 
agents on the sphere is studied, with a constant edge weight 
function for all edges. In particular, in [12], incomplete syn¬ 
chronization is used for accomplishing a flocking behavior, 
where a group of agents moves in a common direction. 
In [13], dynamic agents, which move at constant speed on 
a sphere, are controlled by a state feedback control law 
that steers their velocity vector so as to force the agents 
to attain a collective circular motion; since the agents are 
mass points, the effect of the moment of inertia is not 
studied. In [14], dynamic point mass agents, constrained to 
move on a sphere, are controlled to form patterns on the 
sphere, by constructing attractive and repelling forces; in 
the absence of repelling forces, synchronization is achieved. 
Also, the closed-loop dynamics of these agents are invariant 
to rotations, or symmetry preserving, as those in [11], [12], 
in the sense that two trajectories, whose initial condition 
- composed of position and velocity - differs only on a 
rotation, are the same at each time instant apart from the 
previous rotation. In our framework this property does not 
hold, since our dynamic agents have a moment of inertia, 
unlike the agents in [11], [12], [14], which is another novelty 
of the paper in hand. 

We propose a distributed control strategy for synchro¬ 
nization of elements in the unit sphere domain. The con¬ 
trollers for accomplishing synchronization are constructed 
as functions of distance functions (or reshaping functions as 
denoted in [20]), and, in order to exploit results from graph 
theory, we impose a condition on those distance functions 
that will restrict them to be invariant to rotations of their 
arguments. As a consequence, the proposed controllers can 
be implemented by each agent without the need of a common 
orientation frame. We restrict the proposed controllers to be 
continuous, which means that a synchronized network of 
agents cannot be a global equilibrium configuration, since 
S 2 is a non-contractible set [21], Our main contributions 
lie in proposing for the first time a controller that does not 
require full torque when performing synchronization along 
a principal axis, but rather torque orthogonal to that axis; in 
finding conditions on the distance functions that guarantee 
that a synchronized network is locally asymptotically stable 
for arbitrary connected network graphs, and that guarantee 
that a synchronized network is achieved for almost all initial 
conditions in a tree graph; in providing explicit domains 
of attraction for the network to converge to a synchronized 
network; and in characterizing the equilibria configurations 
for some general, yet specific, types of network graphs. A 
preliminary version of this work was submitted to the 2015 
IEEE Conference on Decision and Control [22]. With respect 
to this preliminary version, this paper presents significantly 
more details on the derivation of the main theorems and 
provides additional results. In particular, the concept of cone 


has been modified, with a clearer intuitive interpretation; the 
proof for the proposition that supports the result on local 
stability of the synchronized network has been simplified; 
further details on the condition imposed on the distance 
functions are provided; additional examples on possible 
distance functions, and their properties, are presented; and 
supplementary simulations are provided which further illus¬ 
trate the theoretical results. 

The remainder of this paper is structured as follows. In 
Section[III] the problem statement is described; in Section IV 
the proposed solution is presented; in Sections [V] and VI 
convergence to a synchronized network is discussed for tree 
and arbitrary graphs, respectively; and, in Section |VlIJ sim¬ 
ulations are presented that illustrate the theoretical results. 

II. Notation 

0„ £ R n and 1„ £ R" denote the zero column vector and 
the column vector with all components equal to 1, respec¬ 
tively; I„ £ R" denotes a column vector in {(x 1; • • • , x n ) £ 
R" : Xi = ±l,Vi = {l,-- - ,n}}; when the subscript n is 
omitted, the dimension n is assumed to be of appropriate 
size. I„ £ R nx " stands for the identity matrix, and we omit 
its subscript when n = 3. The matrix 6>(-)£R 3x3 isa skew- 
symmetric matrix and it satisfies S (a) b = a x b, for any 
a, b £ R 3 . The map II : {x £ R 3 : x T x = 1} ha R 3x3 , 
defined as II (x) = I — xx T , yields a matrix that represents 
the orthogonal projection operator onto the subspace per¬ 
pendicular to x. We denote the Kronecker product between 
A £ R mX " and B £ R aXt by A <8> B £ R msX,lt . Qj ven 
A lt --- ,A„ £ R mXm , for some n,m £ N, we denote 
A = A 1 ® • • • ® A n £ (direct sum of matrices) 

as the block diagonal matrix with block diagonal entries A x 
to A n . Given a,b £ R", a = ±b Oa = bVa = —b; 
additionally, we say a / 0 and b f 0 have the same 
direction if there exists A € R such that b = Aa. We say 
a function / : f2 x ha f2 2 is of class C", or equivalently 
/ £ 0 2 ), if its first n + 1 derivatives (i.e., / (0) , / (1) , 

..., f in >) exist and are continuous on Q,. Finally, given a 
set TL, we use the notation [H\ for the cardinality of T~L. 

III. Problem Statement 

We consider a group of N agents, indexed by the set 
AT = {1, • • • , N}, operating in the unit sphere domain, i.e., 
in S 2 = {x £ R 3 : x T x = 1}. The agents’ network is 
modeled as an undirected static graph, Q = {AC, £}, with AC 
as the vertices’ set indexed by the team members, and £ as 
the edges’ set. For every pair of agents (i,j) £ ACx (AC\{i|), 
that are aware of and can measure each other’s relative 
attitude, we say that agent j is a neighbor of agent i, and 
vice-versa; also, we denote AC C AC as the neighbor set of 
agent i. 

Each agent i has its own orientation frame (w.r.t. an 
unknown inertial orientation frame), represented by 'R, : £ 
SO( 3). Let the unit vector n, £ S 2 be a direction along 




(a) Two non-synchronized rigid bodies. 


(b) Two synchronized rigid bodies (where 
7 ^-i = 7^-2)- 


(c) Two synchronized rigid bodies (where 
7£i / 7^2). 


Fig. 2: In incomplete synchronization, n ri gid bodies, indexed by i = {1, ■ ,n}, synchronize the unit vectors n, = 70h, 

where n* is fixed in rigid body i. In Figs. 2(a)-2(c) rC 
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(u, ,u 2 and u 3 stand for the canonical basis 



(a) Two non-synchronized rigid bodies. 



(b) Two synchronized rigid bodies. 


Fig. 1: In incomplete synchronization, n rigid bodies, 
indexed by i = {1, ■ • ■ , n }, synchronize the unit vectors 
n , 70fi„, where n is fixed in rigid body i. In 
Figs. 1 (a)| 1 (b) rC = -n 2 = ^[1 1 1] T (Uj,u 2 and u 3 
stand for the canonical basis vectors of R 3 ). 


agent’s i orientation, i.e., 11 = 70 il, where h, £ 5 2 

is a constant unit vector, specified in the agent’s i body 
orientation frame, and known by agent i and its neighbors. 
In this paper, the goal of attitude synchronization is not 
that all agents share the same complete orientation, i.e., that 
TZ ± = ■ ■ ■ = 1Z N , but rather that all agents share the same 
orientation along a specific direction, i.e., that 1 ^ = ■ • • = 
n N = ■ ■ ■ = TZ N f\ N . For example, in a group of N 

satellites that must align their principal axis associated to the 
smallest moment of inertia, it follows that, for each i £ A f, 
hi £ S 2 : 3Xt > 0 : 011, = n, with as the satellite’s 
moment of inertia and with A, as the smallest eigenvalue of 
Jp, and that the desired synchronized network of satellites 
satisfies 7Z 1 h 1 = ■ ■ ■ = lZ N h N . Figures [I] and [2] illustrate 
the concept of incomplete synchronization for two agents. 
Notice that agent i is not aware of n, , since this is specified 
w.r.t. an unknown inertial orientation frame; instead, agent i 
is aware of its own direction 11 , - fixed in its own orientation 
frame - and the projection of its neighbors directions onto 
its own orientation frame. 


T ; : R >0 1 —^ R 3 . The rotation matrix 70 : R >0 1 —>■ 50(3) 
evolves according to 

70 (i) = f* (70 (*),«<(*)), (III.l) 

where L : 50(3) x R 3 >->• R 3 * 3 is defined as 

t*(n,u) = ns{u), (IH.2) 

while each unit vector rq : R >0 1 —>■ S 2 evolves according to 
n ; (f) = f„. where f„. : R >0 x5 2 xR 3 £R 3 

is defined as 

f n . (f, n, u) = 5 {Tl, (i)w) n. (III.3) 

The previous results follows from the fact that n, (-) = 
70(-)n, for some constant n ; £ S 2 , and therefore 72 i (f)n i = 
S (1Z i (t)<jJ i (t))'R i (t)h i => n ; (f) = 5 (70 (£)<*>*(£)) n ; (f). 
Finally, the body-framed angular velocity : R >0 1 —> R 3 
evolves according to the dynamics 

4 (70 (*)./<«,(*)) = 70(i)T i (i) (HI-4) 

at 

Ou>i{t ) = J- 1 (-5 J^t) + T,(f)), (III.5) 

and therefore w i (f) = f„. T £ (f)), where : R 3 x 

R 3 i-)- R 3 is defined as 

(«, T) = J- 1 (-5 (w) + T). (III.6) 


Definition 1 : Two unit vectors (n.! , n 2 ) £ (5 2 ) 2 are di¬ 
ametrically opposed if n( n 2 = —1, and synchronized if 
nfn 2 = 1. A group of unit vectors (n l5 • • • ,n N ) £ (5 2 ) N 
is synchronized if nfri, = 1 for alii, j £ {1, • • • ,7V}. 

Problem 1: Given a group of rotation matrices 
(72-j, • • • ,700 : R >0 1 —^ 50( 3) N , with angular velocities 
(wj,-" ,cj n ) : R >0 1 —^ R 3 and moments of inertia 
./,,••• , J N satisfying and design distributed 

control laws for the torques {T, : R >0 1 —> R 3 } ieA r, 
in the absence of a common inertial orientation 
frame, that guarantee that the group of unit vectors 
(nj,-- - ,n„) : R >0 >->■ (5 2 ) N is asymptotically 
synchronized. 


Consider then any agent i £ Af, with rotation matrix 70 : 
R >0 i-> 50(3), unit vector n, : R >0 >->• 5 2 where n ; (-) = 
70(-)rF, body-framed angular velocity u>i : R >0 1 —>■ R 3 , 
moment of inertia ./, £ R 3X3 (,/, > 0), and body frame torque 


For the purposes of analysis, we consider the state x 
(n, uj) := ((n 1; • ■ ■ , n„), (u u ■ ■ ■ , u N )) : R> 0 (->• (5 2 ) n 
(R 3 ) N , and the control input T := (T 3 , ••• , T N ) 
R >0 1 —^ (R 3 ) w ; where x(-) evolves according to x(f) 











f x (f,x(f),T(f)) where 

f x (i,x,T) = (f„(t,n,w),f„(w,T)) xR 3N ,(III.7) 

with f n (f,n,u>) = (f ni (f, n^oti), • • • ,f nj ,(t,n n ,Wj,)) £ 
(R 3 ) N and f„(w,T) = (f^ ( Wl) T\), • • • , f„„ (w W) T w )) £ 
(R 3 ) N . 


Remark 1: Given our choice of state, the corresponding 
vector field in ( [ill] ) depends explicitly on time. The explicit 
time dependency encapsulates the effect of the space orthog¬ 
onal to !!;(•) on the kinematics of h(■). Therefore, in the 
analysis we would have to consider the effect of initializing 
the system at different time instants t 0 £ R, corresponding to 
different initializations of the space orthogonal to n }t Q ). In 
general, that initialization does have an effect on a trajectory 
x(-) of x(f) = f x (i,x(£),T(f)), but, as shall be verified 
later, it does not affect the attainment of the goal defined in 
Problem [T| Therefore, we assume that the initial time instant 
is always 0, i.e., t 0 = 0. 


IV. Proposed Solution 
A. Preliminaries 

We first present some definitions and results from graph 
theory that are used in later sections [23]. A graph Q = 
{A/ - , 8} is said to be connected if there exists a path between 
any two vertices in Af. Q is a tree if it is connected and 
it contains no cycles. An orientation on the graph Q is 
the assignment of a direction to each edge (i,j) £ £, 
where each edge vertex is either the tail or the head of 
the edge. For brevity, we denote N = \Af\, M = \£\ 
and A4 = {1, - - - ,M}. Additionally, and for notational 
convenience in the analysis that follows, consider the sets 
£ = {(i,j) £ Af x Af : j £ Af,}, i.e., the set of edges 
of the graph Q\ and £ = {( i,j ) £ £ : j > *}. For 
undirected network graphs, we can construct an injective 
function R : £ A4 from which it is possible to construct 
a second, now surjective, function k : £ A4 , which 
satisfies n(i,j) = R(i,j ) when j > i and n(i,j) = R(j,i ) 
when j < i. As such, by construction, for every (i, j) £ £, 
n[i,j) = K (ji *)■ since we consider undirected graphs. The 
function «(■, •) thus assigns an edge index to every unordered 
pair of neighbors. The incidence matrix B £ R NxM of Q is 
such that, for every k £ A4 and for (i, j) = R~ 1 {k), B ik = 1, 
B jk = —1 and B lk = 0 for all l £ Af\{i, j}. Finally, for each 
edge k £ A4 and (i,j) = R~ 1 (k), we denote fc n := rp and 
s n := n 3 , i.e., we identify an agent by its node index but 
also by its edges’ indexes (,.n if n is the tail of edge k, and 
s n if n, is the head of edge k ). 

Proposition 2: If Q is a tree, then B T B and (f?®I) T (l?(g)I) 
are positive definite [24], 

Proposition 3: If Q is connected but not a tree, then the null 
space of the incidence matrix, i.e., Af(B), is non-empty, and 
it corresponds to the cycle space of Q (Lemma 3.2 in [25]). 


Denote by C C { 1 , • ■ ■ ,M} the set of indices correspond¬ 
ing the edges that form a cycle. Consider a network with 
n £ N cycles, {C} i={ i, A cycle C, is said to be 
independent if C. L fl C :J = 0 for all j £ {1, • • • ,n}\{«}. In 
Fig. |3(a) a graph with two independent cycles is presented. 
Additionally, we say two cycles C, and C 2 share only one 
edge when | Ci (T C 2 \ = 1 and CV U C 2 contains edges from 
only the following three cycles (in {C} i={1 ... „ } ): C lt C 2 
and C 3 = C, U C.MC, fl C 2 }, with |C 3 | = l^f + |C 2 | - 2. 
Figures 3(b)| and |8(c)| present graphs with two cycles that 
share only one edge. 


Proposition 4: Consider a graph Q with m independent 
cycles, {Ci} i={1 ,... im} . Then the null space of B is given 

by Af{B) = {e = ,e M ) £ R M : e k = ±e i; Vfc ,l £ 

Ci, i = {1, • • • , m}}; and the null space of B g) I„ is given 
by Af(B 0 I„) = {e = (ej,'- ,e M ) £ (K") M : e k = 
±e h \/k,l £ Ci, i = {1, ■ ■ • ,m}}. 

Notice that, for an incidence matrix B £ R NxM , there are 
M edges and these belong to R. On the other hand, for 
the incidence matrix B (g> I„ (with B £ M W><M ), there are 
M edges, but, since the agents operate in an n-dimensional 
space, those edges belong to R". With that in mind, under the 
conditions of Proposition [4] Af(B) = {e = (ei, • • • , e M ) £ 
R" : e k = ±e,,Vfe,Z £ Ci,i = {l,-- - ,m}} means 
that Af(B) is the space where all edges of an independent 
cycle have the same absolute value. On the other hand, 
AT{B (g> I„) = {e = (e 1; • ■ • , e M ) £ (R 71 )" : e k = 
±ei,\/k,l £ Ci,i = {l,--- ,m}} means that Af(B g> I„) 
is the space where all edges of an independent cycle have 
the same direction and norm (or are all zero). 


Proposition 5: Consider a graph Q with n, independent 
cycles, {C} i={1 ... and n 2 pairs of cycles that share only 
one edge, {(C}Cf)} i={li ... ,„ 2} . Then the null space of B 
is given by A f{B) = {e = ,e M ) £ R M : e k = 

±e„Vk,l £ C t ,i = {1, — jnJ} U {e = (e u --- ,e M ) £ 

R" : e fc = ±e„VM G C]\{Cl 0 C},e p = ±e q ,Vp,q £ 
C}\{C] C\C}, i = { 1, • • • , n 2 }}; and the null space of B® I„ 
is given by Af(B®l n ) = {e = (e,, • • • ,e M ) £ (R") M : e k = 

±e„Vfc, l £ Ci,i = {1, ••• U {e = (e ir - - ,e M ) £ 

(R") M : e t = ±e„VM £ Cl\{C]CC*}, e p = ±e q ,Vp,q £ 
C?\{ QnC*},i = {!,-■■ ,n 2 }}. 


The results of Proposition [5] can be interpreted as follows: 
Af(B) is the space where, for each cycle, all its edges 
(£ R), except the one that is shared, have the same absolute 
value; while Af(B 0 I n ) is the space where, for each 
cycle, all its edges (£ R"), except the one that is shared, 
have the same direction and norm (or are all zero). Some 
examples that illustrate Propositions [4] and [5] and details on 
equivalent incidence matrices are found in [Xl] Also, proofs 
of Propositions [4] and [5] are found in Appendix XII These 
Propositions are useful in a later section, where we prove 
that for network graphs that satisfy the conditions of either 
Proposition, the agents converge to a configuration where all 
unit vectors belong to a common plane. 













*(1» 3) = 3 v - / «(l, 5) = 5 

(a) Graphs two independent cycles. 



(b) Graph with one independent cycle and two 
cycles that share only one edge. 


Fig. 3: Two graphs, one with two independent cycles and 
another with one independent cycle and two cycles that 
share only one edge. 



Fig. 4: Three unit vectors n,, n 2 and n 3 in a 30°-cone 
associated to the unit vector n*. 

We now present a definition and some results that will prove 
useful in a later section. 

Definition 2: We say that a group of unit vectors n = 
(n 15 ■ • • , n w ) £ ( S 2 ) N belongs to an open (closed) a £ [0,7r] 
cone, denoted by n £ C(a) (n £ C(a)), if there exists a unit 
vector n* £ S 2 such that n* T n i > cos(ct) (n* T n i > cos(a)) 
for all i £ Af. 

The concept of open o-cone is exemplified In Fig. [4] with 
three unit vectors rij, n 2 and n 3 contained in a 30°-cone 
formed by a unit vector n*. In fact, any group of unit vectors 
contained in the sphere surface region marked in bold is 
contained in a 30°-cone associated to the unit vector n*. 

The following propositions follow from Definition [2] and 
Proposition [27] in Appendix 1 1X | 

Proposition 6: If n = (n 1; • • • ,n N ) £ C(a), for some a £ 
[0, |], then max (j j]eJV 2 (l — nfn,) < 1 — cos(2a). 

This proposition follows immediately after Proposition [27] 

Proposition 7: If, given n = (rij,--- , n N ) £ ( S 2 ) N , 
max (i i)ejV -2 (1 — nJ'nj) < 1 — cos(|a) holds for some 
a £ [0,7r], then n £ C(a). 

A proof of Proposition [7] is found in Appendix [IX] 


B. Distance in S 2 


Consider an arbitrary distance function between unit vec¬ 
tors d : S 2 xS 2 —► Kq", satisfying d( n 2 ) = 0 •<=>■ rij = n 2 . 
We say that d £ V, if, for all (n, , n 2 ) £ 5 2 x S 2 , it satisfies 
the partial differential equation 

S(n l )W^+S(n 2 ) d ^l =0 . (w .n 

a rii on 2 

Let us motivate the introduction of ( |IV-B| ). Given d £ V, 
it follows that, along along a trajectory x(-) of ([III]), (for 
brevity, we omit the time dependencies below) 

d( ni ,n 2 ) = i MM\ i( i, ni ^ i) + M^ r fn2 ( f , n2 , tj2) 

SJ uini s ( ni ) ad( "^; n3) + ^ir 2 s (n 2 ) ad( Z: n " ] 


[Tv-bI 


0 


0 


>l)$( nj 0 ^. (IV.2) 


where we identify an incidence matrix [1 — 1] T corresponding 
to an edge between unit vectors and n 2 . Motivated 
by ( |IV-B| ), let us define two functions. Firstly, consider M 
distance functions d k £ V, one for each edge k £ Ai. Define 
then e k : S 2 x S 2 i->- K 3 as 

dd k (n 3 , n 2 ) 


e fc (n 1 , n 2 ) = S (nj 


Ori! 


(IV.3) 


to be the error of edge k, and for each k £ A4,. Define also 

e : (S 2 ) N ]R 3JV as 

e ( n ) = [ e f (i n , i n ) ••• e^( M n, Al n)] r . (IV.4) 


Define, also, D : ( S 2 ) N >->■ K >0 as 

E k — M 

d k ( k n,- k n), (IV.5) 

fc = l 

named, hereafter, total distance function in the network of 
unit vectors. Note that D( n) = 0 <£=> 3n* e5 2 :n = (1 N (8) 
n*), which means Problem [T| is solved, if along a trajectory 
x(-) of ( fill] ), lirri,^^ D(n(t)) = 0. It follows from ( |IV-B| ) 
that, if d k £ V for all k £ Ai, then, along a trajectory x(-) 
of (gnj), 

D(n(t )) = uj(t)TZ T (t)(B <8 I)e(n(f)). (IV.6) 


Notice that ( | 1V-13 [ ) follows from ( |l V- B [ >, which is the reason 
for imposing such condition. Now, the point in question is 
what types of distance functions satisfy ( |lV-Bj . In Propo¬ 
sition 28 in Appendix, we prove that a distance function 


d(-,-) satisfies dlV-Bl ), i.e. d £ V, if and only if there 
exits / £ C 3 ((0,2),R >0 ), with lim M0 /(s) = 0, such 
that d( n 1 ,n 2 ) = /(I — nfn 2 ). Physically, only distance 
functions that are functions of the angle between two unit 
vectors satisfy ( |IV-B| ) (where the angle between n and 
n 2 can be defined as 0(n 1 ,n 2 ) = arccos(nfn 2 ), with 
9 : S 2 x S 2 i —> [0,7r]). Moreover, distance functions that 
satisfy ( |IV-B| > are invariant to rotations of their arguments, 
i.e., d(7£ni, 7 ^ 112 ) = d(ni,n 2 ) for any 1Z £ 50(3). 
This property guarantees that the proposed controllers can 
be implemented without the need of a common inertial 
orientation frame, a property that is verified later. 


To summarize, there are M distance functions d k £ V, 
one for each edge k £ Ad, and, therefore, there are M 
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Fig. 5: Relation between properties of /(•) and the classes 
it belongs to. 


f k £ C 1 ((0, 2),K >0 ), one for each edge k £ Ai, such that 
d fc (n 1; n 2 ) = f k ( 1 — nfn 2 ). For reasons that will become 
apparent later, we restrict the previous functions /*.(•) to 
satisfy some more properties. 


Definition 3: Consider a function / £ C 2 ((0, 2), R >0 ), sat¬ 
isfying i) f(s) > OVs G (0, 2), ii) lim s _ >0 + /(s) = 0, and Hi) 
limsup a ^ 0+ < oo. Denote f 2 := lim s ^ 2 - /(s) 

and fg := lim s ^ 0 + /'(s). We say 


. / G V 0 if /' = 0 and / G V B if /' ? 0, 

• / G V°° if f 2 = oo, and / G if / 2 < oo, 

. / G V° if / £ V & A lim s _ >2 - f{s)s/^~s = 0, 
. / G V_° if / £ V* A lim s _ >2 - f (s)y/2-~h ± 0, 

• / £ T 7 if /(•) is of any of the previous classes. 


Figure [6] illustrates the different classes introduced in Defi¬ 
nition [3]while Fig. [5] illustrates how the properties that /(•) 
satisfies affects the classes it belongs to (see Remark [33] 
in Appendix). In [20], the notion of reshaping function is 
introduced, whose definition is within the same spirit as 
that of Definition [3] For the rest of this manuscript, we 
assume that, for each edge k G M, there exists a function 
d k (n u n 2 ) = / fc (l - nfn 2 ) where f k G P; in particular, 
f k (-) plays the role of an edge weight. In e.g. [11], [12], 
/k(s) = a k s and /'(s) = a k , for all k G M (a k is 
the weight of edge k and it is denoted by a, :l in [11], 
where (i. j) = ic -1 (fc)). As such, ( |IV-B] i and ( |IV-B| > can 
be rewritten as D : 17° i-> R >0 (see ( |IV-B| )) 

E k — M 

f k (l- k n T k n) (IV.7) 

fc = l 

and e fc : {(n 1; n 2 ) G S 2 x S 2 : nfn 2 ^ -1 if f k £ V°) (->■ 
R >0 as 

e^n^n,) = S (nj 9fk( ) h f l 1 " 2) = /((1 - nfn 2 )S (njn^ (IV.8) 
=> ||e t (n 1 ,n 2 )|| = ft{s)^s(2 - s) |, =1 _„r„ 2 , 


respectively. Note that e fc (n 1 ,n 2 ) is well defined for all 
(n^iia) G S 2 x S 2 if f k G V°, and if f k G V°, note that 
lim e fc (n 1 ,n 2 ) = lim f' k {s)^s(2-s) lim n|^gf|y 

ns-tn! s-»2- na-t-nj M \ -II 

does not exist. Note that the total distance function ( |1V-B| ) 
depends on /*.(•), for all k G A4, while ( ]1 V-B[ > depends on 
/'(•), for all k G Ai. As such, a distance function may or 
may not be defined when two unit vectors are diametrically 
opposed, depending on whether / G V°° or / £ V 
similarly, an edge error may or may not be defined when 
two unit vectors are diametrically opposed, depending on 
whether / G V° or / G V°. Thus, the domains of ( |IV-B[ ) 
and of ( ]IV-B| ) depend on the classes / fc (-) belongs to, for all 
k G A4. In particular, D : i-> R >0 and e : fl® i-> R 3JV , 


where 

nf ={n G ( S 2 r ■ + —1 V/ fc G P”}, (IV.9) 

ni ={n G (5 2 ) N : k n T - k n ± -1 V/ fc G V 0 }, 

and where we emphasize that 17“ C 17^, since f k G V°° => 
f k £ V" (see Fig. [ 5 }. These domains play a role later on, 
since D (■) is used in constructing a Lyapunov function, while 
e(-) is used in constructing the control law. As such, the 
Lyapunov function can be well defined, while the control 
law is not, while if the control law is well defined, so is the 
Lyapunov function. Consequently, it is important to guar¬ 
antee that along trajectories of the closed-loop system, the 
control law is well defined. Additionally, notice that ( |IV-B| ) 
provides some insight on why we denote e fc (-,-) as edge 
error of edge k. Indeed, if / fc G 'PjVfc G A7, it follows that 
e k ( k n , s n) = 0 implies that fc n = ±j;.n, i.e., it implies that 
the neighbors that form edge k are either synchronized or 
diametrically opposed. Moreover, if f k G V V/.' G Ai, the 
distance between unit vectors is supremum when two unit 
vectors are diametrically opposed, i.e., for each k G Ai, 
(denote (7 = {(n 1; n 2 ) G S 2 x S 2 : nfn 2 = —1}) 

sup d t (n 1 ,n 2 )= sup / k (l - nfrij) = lim f k (s) =: d™ x . 

( n l’ n 2)£f2 (r»i,n2)GS7 

For convenience, denote 

d nul ' := iriind:: ,a L (IV.10) 

k&M 

which plays an important role in this and the following 
sections. 

Proposition 8: Consider the total edge error in ( | IV- B [ > and 
the total distance function in ( |IV-B| ). Consider Vf : R >() 
f7f as Q,' n (D) = {n G fl° : D(n) < D}, where D.' n (D) is 
compact for all positive D (where refers to a set-valued 
function). Then, it follows that 

\/D < d min , max ||e(n)|| < 00 , 

nel!^(D) 

and that there are no diametrically opposed neighbors, i.e., 
|{g £ A7 : Vn G ,n T ,-n = —1}| = 0. If f k G 

V° for all k G Ai, it follows that max„ en e ||e(n)|| = 
max ne(s 2 ) iv ||e(n)|| < 00 ; moreover, given D < pd min 
for some p G A7, it follows that there are at most p 1 
diametrically opposed neighbors, i.e., |{g G Ai : Vn G 
,n T 9 -n = — 1 }| < p - 1 . 

Proof: Recall Definition[3] ( |l V-B[ ) and ( | 1V- B[ ). It follows 
that e(-) and D(-) grow unbounded only if two neighbors 
are diametrically opposed (however, this is not sufficient). 
Recall also that 17“ C O", meaning that e(-) can grow 
unbounded while D(-) remains bounded. Finally, notice that, 
if n G ( S 2 ) N is such that there are q diametrically opposed 
neighbors, it follows that D(n) > qd"" n . Consider then a 
D < d""” < 00 . Then, for all n G it follows that 

D( n) < D =► Vfc G Ai, fc n T j.n > 1 - (D) > -1, 

which means that two neighbors are not arbitrarily close to 
a configuration where they are diametrically opposed and, 
therefore, max nen / (fi) ||e(n)|| is bounded. For the second 
part of the Proposition, consider that f k G V" for all k G Ai. 
Then f7“ = f7f = (S 2 ) N , which means max ne(1 ^ ( n) || e ( n )|| 

































Fig. 6: Four functions belonging to different classes as 
introduced in Definition [3] (from top to bottom in legend) 

f(s) = s, f(s) = it 2 arccos 2 (l — s), 
f(s) = tan 2 (0.5arccos(l — s)) and 
/(s) = 0.25( v / s(2 - s)(s - 1) + arccos(l — s)). 

and max nel) / (S) fl(n) are bounded (since the domain is 
compact, and the functions are continuous). Consider then a 
p £ AA and a D < pd min < oo. Since D( n) > qd min , then, if 
there are q diametrically opposed neighbors in n £ ( S 2 ) N , it 
follows that |{q € AA : Vn £ f,n T 9 n = — 1}| < p— 1. 


C. Solution to Problem [7] 


In this section, we present the controllers for the torques 
of each agent i. For each agent i, we design a controller that 
is a function of \Af r + 1 measurements: \Af r measurements 
corresponding to the distance measurements between agent i 
and its \Af r neighbors, and 1 measurement corresponding 
to the body frame angular velocity. More specifically, we 
assume each agent i measures TZJ (f)n d (t) = 7 ZJ (f)7?. j (f)n J , 
at each time instant t > 0, and for each j £ Af ; physically, 
this means that agent i knows n ; (the constant unit vector 
that it is required to synchronize with), and that it can 
measure the projection of this unit vector on its orientation 
frame; each agent i must also measure which does 

not require an inertial reference frame. For convenience 
denote A/” = {*!,••■ ,}, and, given n, £ S 2 , denote 

,n I|Aft| ) <E S 2|JVil : nJn H ± -1,VZ £ 
{1, • • - , |A/71} A £ V 0 } which provides the domain 

where the control law for agent i is well defined (recall that 
if / fc £ V°, for some k £ AA, then ( |IV-B| ) is not defined 
when two unit vectors are diametrically opposed). We then 
propose, for each agent i £ Af, the following decentralized 
control law Tf : = (( IZJn i± , 

f xR 3 gK 3 




Tf = -er^) - ( n i,^n i( )jV.ll) 

with a £ £ (see Definition [5] in Appendix). The timed 
control laws for each agent i £ Af are then T ; : R >0 i—>• R 3 
as 


T i (f) = Tf((^(f)n, 1 (f) I ... ,^(t)n iM| (t)) J w 1 (t))OV.12) 

The proposed torque control law exhibits the following 
properties. The controller function in ( |IV-C| ) is decentralized 
in the sense that it does not depend on the measurement 
of the global state x(-). Also, ( |IV-C[ i can be implemented 
without the knowledge of an inertial reference, since mea¬ 
suring yi(t) := HT (i)lZi (t)hi at every time instant t > 0 
and for all Z £ {1, • • • , |Wj |} only requires the measurement 
of the projection of n,, in agent’s i body orientation frame; 
while y' A ..| +1 (t) := OJ^t) is also measured in agent’s i body 
orientation frame. Finally, notice that ||T i (-)|| < cr m “ + 
\Afi\ max jeAf . sup 0<a<2 /' (iJ) (s). As such, the proposed con¬ 
trol law, for each agent i, can be implemented with bounded 
actuation provided that cr max < oo and that £ V° for 


all j £ Afi. 


Notice that IM‘> Oil < IM‘> Oil < 

\/M||e(-)||, and, therefore, for any 7?.; £ SO( 3), and for 

all x G x R 3JV , 


|T° i ((72.fn il , ■ ■ ■ ,WTn i|Vj| ),w i )|| < 0-j|u>4 + v / M||e(n)PV.13) 

a relation to be used later. We can now present the complete 
control law T c! : R >0 x f2® x R 3JV i-x R 3JV , 


- T 1 ((7Z[ (fjnjj,■ ■ ■ ,7e> 1|Ari| ), Wl ) 

T c ‘(t,x)= : 

_T N ((TZ'^(t)n Nl , • • • ,7^.^n JV| ^ | ),a; i v)_ 

= -[cr T (uj 1 ) ■ ■ ■ <t t {u) n )] t - U T (t)(B ® I)e(n), (IV. 14) 


where TZ(-) = 7£i(-) ® • • • ® 7^ JV (-) (see Notation). For the 
remainder of this paper, we dedicate efforts in studying the 
equilibria configurations induced by this control law (for 
different types of graphs), their stability, and what is the 
effect of the chosen distance functions. Notice that ( |IV-C| ) is 
defined on R >0 x xR 3JV . As such, when f k £ V°\/k £ A4, 
= ( S 2 ) N , and the analysis is simplified; when, however, 
3 k £ Ad : /i £ ?*, C ( S 2 ) N (where D® is open), 
and it is necessary to guarantee that a trajectory x(-) never 
approaches the boundary of 


D. Lyapunov Function 


In addition to the total distance function of the net¬ 
work ( |IV-B| ), let us also define the total rotational kinetic 
energy of the network, as H : R 3N i—>• R >0 , where 

1 v- yi=N 

H{u) = 2E h u > 3 (IV15) 

and which satisfies f x (-,x, T) = f u (u>,T) = 

E::f = uj t T, for all (w,T) £ R 6JV . 

Combining ( |IV-B[ ) and ( |l Y’-D[ i, consider then the Lyapunov 
function V : fl^ x R 3JV hg R >0 , 

V(x) = D(n) + H(u), (IV. 16) 

and the function W : R 3JV i-x R >0 defined as 

fx(‘j x ,T c! (-,x)) 

^ - u> T K T {-)(B ® I)e(n) - w r T ci (-, x) 

(IV. 17) 

and it follows that, along a trajectory x(-) of x(t) = 
f x (7,x(f),T ci (f,x)), V(x(t)) = -W(v(t)) < 0,Vt > 0. 
Note that, for every i £ Af, 

IM^Tf(-))||? x ^(Tf(.) + A m „(J,)||o;,|| 2 ) 

^ A ( cr 'H a2 -ll + ^ll e ( n )ll + '’w(T) ll^.ll 2 ) =: /”(x).(IV.18) 
























Moreover, along a trajectory x(-), (for brevity, below, we 
omit the time dependency of the state) 

• dW (x) T 

= ^ fx(f) x, T ol (t, x)) 


= S 1=1 


duJt 




for which it follows that Vi > 0 


|W / (w(i))| < (°A + a ') T l|Wi(i)||/”(x(i)).(IV.19) 

z — J \ =i 1 

Along the same trajectory, it also follows that, for every i G 
TV (we again omit the time dependency) 


u), = — J, 1 (5 (tjj) JjO;, + S (u»i) JjOij + _D<t(oj,)<A) + 

J-'ftfSfn,) £ (/'(.)I -/"(>X) 5 (n^^-^OV^O) 

It follows from dTV-Pl ), ( flTODj i and ( flTODj i that if 

su Pt>o ||e(n(i))|| < 00 and sup t > 0 ||u> i (i)|| < oo, then 
sup t>0 \W(u>(t))\ < oo and sup t>0 ||d> i (i)|| < oo; this in 
turn implies that W (<*>(•)) and are uniformly continu¬ 
ous, which plays a role in proving that lirri,^, W= 0 
and that lim^^ ujft) = 0. 


Proposition 9: Consider the vector field (JTTTji, the con¬ 
trol law dlV-C), and a trajectory x(-) of x(i) = 
f x (i,x(i),T ci (i,x(i))). If x(0) G = {x £ n° x R 3JV : 
V(x) < d min }, then, along the trajectory x(-), lim t _ >00 (f? ® 
I)e(n(f)) = 0 and lini,^^ oj(i) = 0. 


We refer to the proof of Proposition [TO] which provides a 
result for a more general control law than that in ( |IV-C| >, and 
which is described in the next section. 


E. Constrained Torque 

A natural constraint in a physical system is to require 
the torque provided by agent i to be orthogonal to n,. In 
satellites, thrusters that provide torque along n ; might be 
unavailable; also, controlling the space orthogonal to n, 
can be left as an additional degree of freedom, in order 
to accomplish some other control objectives. However, the 
control laws proposed in ( |IV-C| ) require full torque actuation, 
in particular, ( |IV-C| > requires each agent to provide torque on 
the plane orthogonal to n,. Indeed, since nfe fc (n 1 ,-) = 0, 
Vn, G S 2 ,\/k G AA (see ( jIV-B| i), it follows that, for all 
i G A f, 

nfTfO,^) = rH<x(u>i), 

which is not necessarily 0. In short, previously, we provided 
control laws Tf : Q n . x R 3 i-> R 3 which require full 
torque by each agent i G A f, and in this section we provide 
constrained control laws T f : Sl n; x R 3 i->- {z G R 3 : 
z T iL = 0}, i.e., control laws which do not require torque 
along hi. Let us anticipate a future result by announcing 
that the constrained control law can only be used by agent 
i G Af when the unit vector to be synchronized by agent 
i G A f, namely bp, is a principal axis of that agent (i.e., when 
n, is an eigenvector of ,/,). Consider then Tf : (i/ 4 , wf) = 
,TZfn iw ),u}i) G Q„, x R 3 {z G R 3 : 


z T n t = 0} defined as 

Tf (ia, u)i) = n(n i )Tf(v„w i ) 


E3 


_ \ ^ I — | A/”~ | 

- cr ( n (n i )w i )e s(M 


(h,,^;'n, i )(.IV.21) 


Additionally, consider a partition of Af, i.e., £u£ = Af with 
£p\£ = 0; where £ is a subset (possibly empty) of the agents 
whose unit vector to synchronize is an eigenvector of their 
moment of inertia, i.e., £ C {i £ Af : 3A i s.t. J,n, = A,n;}- 
Then we propose the complete control law T' ' : (t. x) G 
R >0 x x R 3JV i->- R 3JV defined as 


’ (IV.22) 


I (e, ® l 3 ) T T-(t,x) = TfffKfWn.,, ■ ■ ■ Vi 6 

\(e f ® l 3 ) T T“(i,x) = T?((KJ(t) n.,, • ■ ■ V, „ 

i.e., for agents whose unit vector to synchronize is a prin¬ 
cipal axis, either control law ( |IV-C| ) or ( IV-E) is chosen, 
and, for all other agents, control law ( |IV-C| l is chosen. As 
such, agents whose unit vector to synchronize is a principal 
axis have an option between using full torque control or 
constrained torque control. The disadvantage with the control 
law in ( |IV-E| > is that, along a trajectory of the closed-loop 
system, and for all i G £, lim^^ h'ujff) is not guaranteed 
to exist and be 0; i.e., an agent that opts for ( |IV-E| ) can 
asymptotically spin, with non-zero angular velocity, around 
h t (nonetheless, we can guarantee that sup t>0 ||u> i (f)|| < 
oo => sup t>0 n'w.lT) < oo, i.e., if an agent applies ( |IV-E| ) 
it never spins infinitely fast around its principal axis n,). 


Since ||n(n,)u; i || < ||,Vn^ G 5 2 ,a; i G R 3 , it follows 

that ||f„. Tf (•))!! < /”(x), with /”(•) as defined 
in ( |!V-D| . Similarly to ( fTVd5| , denote W : R 3JV ^ R> 0 
as 

^ (w)=_ ^T fx(-,x,T d (f,x)) 

= J2 u} > (T ( UJi ') + H^ Tn (TO) < 7 ( w i). 

iec jec 

and it follows that, along a trajectory x(-) of x(t) = 
fx(£x(f),T ci (f,x)), H(x(f)) = -W(u>(t)) < OVf > 0. 
Moreover, along a trajectory x(-), (for brevity, below, we 
omit the time dependency of the state) 

= f x (t,x,T e ‘(f,x)) => 

i — N 

=> \W(u{t))\ < +a') T ||w i (f)||/“(x(f)). aV.23) 

1 = 1 

Proposition 10: Consider the vector field ( |f!![ >. the con¬ 
trol law ( |IV-E| >, and a trajectory x(-) of x(t) = 
f x (f, x(f), T oi (f, x(f))). If x(0) G n° x = {x G 0° x 
R 3JV : V(x) < d min }, then lim t ^ 00 (B (g> I)e(n(t)) = 0, 
\irn. t ^ 00 u}i(t) = 0 for i G £ and lim^^ n (n^) ur, (f) = 0 
for j G £. Moreover, sup t>0 |nJWj(f)| < oo for j G £. 

Proof: For brevity, we say f : R >0 i—>• R" is bounded, 
if sup (>0 ||f(f)|| < oo; we say f : R >0 i-> R" converges 
to a constant, if 3f“ G R" : lim t _> 00 f(f) = f“. Let 
us provide a brief summary for the proof. First, we prove 
that, along a trajectory x(-), ||o;(-)|| and ||e(n(-))|| are 
bounded. This, in turn, guarantees uniform continuity of 









































V(x(-)) and of u>(-). And finally, since both V(x(-)) and 
<*>(•) converge to a constant, we invoke Barbalat’s lemma 
(see [26], Lemma 4.2) to conclude that e(n(-)) converges to 
the null space of B (g> I. Recall then the functions in ( ]l V-D[ ) 
and ( |IV-D) . Since V'(x(-)) < —W"(<*>(•)) < 0, it follows that 
V(x(-)) < V(x(0)) < d mi ". Therefore D( n(-)) < d min and 
H(uj(-)) < d min . From D( n(-)) < d min , it follows, with the 
help of Proposition [8j that e(n(-)) is bounded; while from 
if (<*>(■)) < d min , it follows that u>(-) is also bounded. From 
boundedness of e(n(-)) and u>(-), it follows that T ci (-,x(-)) 
is bounded (see ( |IV-E| ) and ( |IV-C| >); that ||w i (-)|| < /”(x(-)) 
is bounded (see ( IV-D[ |); that |V r (x(-))| = |VF(o;(-))| is 
bounded (see ( |IV-E i); and, finally, that w,(-) is bounded 
(see ( |IV-D| >). The previous conclusions imply that V (x(-)) 
and that <!>(■) are both uniformly continuous. Since V(-) > 0 
and V(x(-)) < -W(u>(-)) < 0, it follows that V(x(-)) 
converges to a constant; by Barbalat’s lemma, uniform con¬ 
tinuity of V(x(-)) then implies that V(x(-)) = —W(&(■)) 
converges to 0. As such, it follows from ( |IV-D| >, that a; i (-) 
converges to 0, for all i £ C, while II (n,) w,(-) converges 
to 0, for all j £ £\ also, notice that 


limn^Jw^t) = 0 => ImiK(f) - n^njw^f))) = OfIV.24) 

Let us now study agents in C and £ separately. Also, 
for convenience, and with some abuse of notation, denote 
Tf(£) = (e 4 <g) l 3 ) T T ci (f,x(f)), for i £ C, and T f(t) = 
(e,,-C)l 3 ) T T ci (f, x(f)), for j £ C. For i £ C (for which ( |IV-C| ) 
is the chosen control law), and again by Barbalat’s lemma, 
convergence of <*>;(•) to 0 and uniform continuity of uj,(-) 
imply that u;^-) = f u (<*£;(■), Tf (•)) converges to 0; since 
<*>;(■) converges to 0, so does Tf (-,x(-)) (see [ill])- Now, for 
j £ £ (for which ( |l V-L| i is the chosen control law), and once 
again by Barbalat’s lemma, convergence of II (n j w, (-) to 
0 and uniform continuity of ^(11 (rij) <*>(£)) = II (n^) uj(t) 
implies that II (n,) di(-) converges to 0, and therefore 


liin n(n J )tj J .(t) = 0 => lim(Wj(£) - n^nTw.,(t))) = 0(1 V.25) 

Now, recall ® where J^ jt) = -S {u^t)) + 

T f(t), and, from ( |1 V-C[ i and ( |IV-L[ ), it follows that 

liiii + S (flj) fnfnjcj^t)) 2 - Tf(i)) = 0 

'' i " 3 lim ;n (nju^i)) - T f(t)) = 0 (IV.26) 

If we take the inner product of ( |IV-E| > with n,, and 
since T'v'(-) _L n 2 , it follows that lim t _ >00 (nj’tb i (£)) = 
linit-foo nJ’Tj’(-) = 0. As such, it follows from ( |lV-L[ ) 
that T c '(-) converges to 0 for all j £ £. Now to 
summarize, recall that, for all i £ C, both <*>;(•) and 
Tf(-,x(-)) converge to 0, which implies, from ( |IV-C| ), 
that (n 4 , lZ i (-) T n il (•)) converges to 0. On 

the other hand, for all j £ £, both II (n ;l ) oj :) (•) and 
T c! (-,x(-)) converge to zero, which implies, from ( ]lV-L) , 
that (n.-^(')X(’)) converges to 0. All 

together, it implies that (B (g> I)e(n(-)) converges to 0. 
Finally, sup t>0 |n T u; j (£)| < oo since <*>(■) is bounded 

< cP"). ■ 


Notice that d""", defined in (|IV-B[i, is a design parameter, and 


therefore, the domain of attraction in Proposition p~0] can be 
made larger by increasing this parameter. More specifically, 
d"" u increases the domain of attraction in the state space 
related to u>, which is clearer in the next corollaries. 


Corollary 11: Proposition |9j holds if ; 

< 1 — r. 


D( n(0)) 

d mm 


H{u>( 0 )) 
d mm 


< 1 and 


Proof: If the Corollary conditions hold, then V (x(0)) = 
D(n(0)) + H(uj(0)) < d mi " - H{u( 0)) + H(u(0)) = d min . 


Corollary m states that if the total kinetic energy is small 
w.r.t. d""", and if the total distance in the network is small 
w.r.t. 1 — r, asymptotic synchronization is guaranteed. 

Corollary 12: Proposition [ 9 ] holds if r := < 1 and 

if 


d fc ( fc n(0), fcii(O)) 1 — r WI _ 

y v , y " < —— ,Mk £ M. (IV.27) 

cH n M 

Proof: If ( fl2| holds, then H(n(0)) = 

E fce £ d fc(fc n (°)»i n ( 0 )) < Mmax teA1 d k ( k n(0), sn( 0 )) < 
rfm' n [I — r), in which case. Corollary 11 holds. ■ 


Corollary 12 states that if the total kinetic energy is small, 
and if all neighbors are close, then asymptotic synchroniza¬ 
tion is guaranteed. 


Corollary 13: Proposition [ 9 ] holds if r := < 1 and 

if 

n ( 0 ) e C ^5 arccos ^1 - min f- 1 , (IV.28) 

where d min = min fcg7M lim,^ 2 f k (s). If f k £ V°° for all k £ 
A4, then d min = 00 and ( p~3| > reduces to n(0) £ C(|). 


Proof: By invoking Proposition [6j it follows that ( fT3| ) 
implies that d k ( k n(0), s n(0)) < A(min/- 1 (d min i=P)) < 

l GM 

d""" 'ffj- for all k £ M, and therefore ( fl2] ) holds. ■ 


Corollary [13] states that if the total kinetic energy is small, 
and if all neighbors are initially contained in a small cone, 
then synchronization is guaranteed. Moreover, if d'""‘ = 00 
and if all neighbors are initially contained in an open j -cone, 
then synchronization is also guaranteed. 


Proposition 14: Consider the vector field ( |fII[ >, the con¬ 
trol law (|IV-E| i, and a trajectory x(-) of x(t) = 
f x (£,x(£),T ci (£,x(£))). If f k £ V° for all k £ M, then 
for all x(0) G (S 2 ) N x K 3JV , lim^^^ 0 I)e(n(f)) = 0, 
lim^^ u>i(f) = 0 for i £ C and lim^^ II (n 2 ) ax, (f) = 0 
for j £ £; additionally, if x(0) G f2° = {x G ( S 2 ) N x R 3W : 
3p £ M,V (x) < pd min }, then no more than p— 1 neighbors 
are ever diametrically opposed, i.e., sup t>0 \{q £ A4 : 
,n T (f),-n(f) = -1}| <p- 1. 


Proof: Notice that if f k £ V° for all k £ M then 
= = (‘^ 2 ) JV > which is a compact set. Since e(-) 

is continuous in 12®, it follows that ma,x nfQ c ||e(n)|| < 
00, and, therefore, ||e(n(-))|| is bounded regardless of the 
trajectory x(-). To conclude that lim t _ >00 (f3(g)I)e(n(f)) = 0, 
lim^oo (jjft) = 0 for i £ C and lim^^ II (n 2 ) ox, (f) = 0 
















































for j £ C, it suffices to follow the same steps as in the proof 
of Proposition [TO] For the final statement in the Proposition, 
consider x(0) £ f2“ = {x £ f2 x : 3p £ A4,V(x.) < 
pd min }. Since, along a trajectory x(-), Z2(n(-)) < V(x(-)) < 
V(x(0)) < pd"'"\ it suffices to invoke Proposition [8] with 
D = V(x(0)), and the Proposition’s conclusion follows. ■ 

Denote x) := f x (f, x, T cl (t, x)) as the closed-loop 
vector field. Note then that QQ = {x £ ( S 2 ) N x R 3JV : 
Vt > 0. if (t, x) = 0} provides the set of all equilibrium 
points, and moreover {x £ (5 2 ) JV x R 3JV : (B <8> I)e(n) = 
0, Wi = 0 for i £ C, II (n./) u>j = 0 for j £ £} C (2(). As 
such. Propositions [TO] and [l4| imply that, under the respective 
Propositions’ conditions, a trajectory x(-) converges to the 
set of equilibrium points. Note also that [(l w ®n* ) T -] T £ 12^ 
for all n* £ 5 2 , i.e., all configurations where all agents 
are synchronized are equilibrium configurations (agents are 
synchronized and not moving, or agents are synchronized 
and spinning around their principal axis). Finally, notice that 
since e(Sn) = e(n) for all S £ {1^ ® 7 Z £ R 3NX3JV ; 1Z £ 
50(3)} and for all n £ f2®, it follows that has geometric 
isomerism [14]; i.e. [n T -] £ 12^ => [Sn T -j £ which 
means that for every equilibrium configuration, there exits 
infinite other equilibria configurations which are the same 
up to a rotation. In Section [V] for tree graphs, we show that 
fis composed of configurations where agents are either 
synchronized or diametrically opposed; while in Section [VT] 
for graphs discussed in Propositions [4] and [5] we show that 
12^ is composed of configurations where agents belong to a 
common plane. In light of these comments, it follows that 
Corollaries [TTpTI provide conditions for when a trajectory is 
guaranteed to converge to a configuration where all agents 
are synchronized, and not any other configuration in 1 P ; 
in particular, if the initial kinetic energy is too large with 
respect to d'""\ the agents may escape to other equilibria 
configurations other than synchronized ones. 

Remark 15: Consider the closed-loop vector field 
f°'(f,x) = f x (t,x, T c! (f,x)). Additionally, consider the 
alternative state x(t) = diag(I 3N , TZ(t))x.(t) 44 x(t) = 
diag(I 3JV , 1Z T (t))5t(t), which evolves according to 

x(i) =diag(0,7t(f))x(i)+ 

diag(I 3N 1 x(t)=diag(I 3 AT,7i :r (t))x(t) 

= diag(I 3JV ,7?.(f))f x ‘(f, diag(I 3W , 7^ :^ (f))x(^)) 
=:F c! (f,x(f)) 

The difference between x(-) and x(-) is that all quantities 
in x(-) are expressed in the inertial reference frame, while 
in x(-) the angular velocities are expressed in the body 
reference frame of each agent. If = jj for all i £ Af, 
then for any S £ {TZ ® l 2N £ R 6Nx6JV : TZ £ 50(3)}, it 
can be verified that f}'(-,Sx) = Sf5(-,x), which implies 
that the closed-loop dynamics of the agents are invariant to 
rotations; i.e., given two initial conditions x^O) and x 2 (0) 
satisfying x 1 (0) = Sx 2 (0), it follows that x 3 (f) = Sx 2 (t) 
for all t > 0. This is the case, because when J t = j,\, 
then reduces to a second order integrator. A similar 


result has been reported in [14], where the agents are mass 
points (i.e., agents without moment of inertia). However, in 
our framework, where in general ./, / j J for some i £ TV", 
invariance of the closed-loop dynamics to rotations does not 
hold due to the term 5 (u>,) d,u3, in ([ill]). 


V. Tree Graphs 


Let us focus first on static tree graphs. For these graphs. 
Proposition [2] states that A’(/iQl) = {0}. In this section, we 
quantify the domain of attraction for synchronization to be 
asymptotically reached, i.e., we construct a domain 12“ such 
that if x(0) £ 12“, then all trajectories of ( |I 11 [ > under control 
law ( |IV-E| > asymptotically converge to a configuration where 
all unit vectors are synchronized. Later, we construct another 
12 “, for graphs other than tree graphs, which is smaller in 
size, and we quantify how much smaller it is. 

Theorem 16: Consider a static tree graph, the vector 
field ( [ITT] ), the control law ( |IV-E| ), and a trajectory x( ) of 
x(f) = f x (i,x(f),T ci (f,x(i))). If x(0) £ 12“ = {x £ 
12° x R 3N : V(x) < d min } then synchronization is asymp¬ 
totically reached, i.e., lim^^ (rq (t) — n ; (f)) = 0, for all 
(i,j) £ Af 2 . If f k £ V°° for all k £ A4, then d min = 00 
and synchronization is asymptotically reached for almost all 
initial conditions in (5 2 ) N x R 3JV . 


Proof: Under the Theorem’s conditions, we can invoke 
Proposition 10 to conclude that lim t _ +00 (.B ® I)e(n(f)) = 


0 and Proposition [8] to conclude that two neighbors are 
never arbitrarily close to a configuration where they are 
diametrically opposed. Since Af{B (g> I) = {0}, it follows 
that lim^^i? ® I)e(n(f)) => linq^^ e(n(f)) = 0. As 
such, and since two neighbors are never arbitrarily close 
to a configuration where they are diametrically opposed, it 
follows that all unit vectors converge to one another. For 
the second part of the Theorem, notice that, if d'""' = 00, 
then f2“ = f2° x R 3iV . Since (2° x R 3JV \{(5 2 ) JV x R 3JV } = 
{x £ (5 2 ) N x R 3JV : V/ fc £ 'P 00 , fc n T s n = -1} is a set 
of zero measure in the space of all initial conditions, i.e. 
(5 2 ) N xR 3lv , synchronization for almost all initial conditions 
is guaranteed when d min = 00. ■ 


Notice that in Theorem 16 increasing <2 mi " enlarges the 
region of stability, and it yields the almost global stability 
result for d min = 00. However, a similar result for other 
graphs, than tree graphs, is not presented in this manuscript. 


Example 1: Consider the distance functions d(ni,ri 2 ) = 
/(1 — n}n 2 ) where f(s) = a ^ arcco3 ( 1 ~ 5 ) ^ > with a > 0 
and a > 2. For these, d max = a, f £ V * and / £ V°\ 
also / £ V 0 . Suppose f k (s) = f(s) for all k £ Ad, and 
for some a and a. Invoking Corollary [13] it follows that if 
r := < l and 


n(0)£C - 


1 — r 
M 














then Theorem’s 16 conclusions follow. Notice that by in¬ 
creasing a convergence for arbitrary initial values of rota¬ 
tional kinetic energy can be guaranteed; on the other hand, 
by increasing a we can increase the size of the cone where 
the agents need to initially be contained (up to a cone C (f ))■ 
We emphasize that the domain of attraction in Theorem [l6| 
is however larger, in the sense that there are more initial 
conditions which do not satisfy the previous conditions, but 
for which synchronization is still guaranteed. 


Example 2: Consider the distance function d(n 1 ,n 2 ) = 
/(l — n^rn) where /(s) = 2 a tan 2 (0.5arccos(l — s)), for 
which d max = 00 and / £ V°° (and therefore / £ V°). 
If f k (s) = f(s) for all k £ A4, then convergence to a 
synchronized network in a tree graph is guaranteed for almost 
all initial conditions (see Theorem m- 

Theorem 17: Consider a static tree graph, the vector 
fields ( |11I| ), the control law ( |I V-L£[ ). and a trajectory x(-) of 
x(t) = f x (i, x(i),T ci (f, x(f))). If f k € V° for all k £ M, 
and x(0) £ = {x £ xR M : 3p £ M, V(x) < pd min } 

then the group of unit vectors converges to a configuration 
where no more than p — 1 neighboring unit vectors are 
diametrically opposed. 

Proof: Under the Theorem’s conditions. Proposition [l4| 
can be invoked. Additionally, since Af(B (g) I) = {0} in 
a tree graph, it follows that lim t (g> I)e(n(t)) => 

lim^^ e(n(t)) = 0, which implies that all neighbors are 
either synchronized or diametrically opposed. Since, by 
Proposition [14] there are at most p— 1 diametrically opposed 
neighboring unit vectors, it follows that the group of unit 
vectors converges to a configuration where no more than 
p — 1 neighboring unit vectors are diametrically opposed. ■ 


shown to be unstable. For proving this statement, it suffices 
to find initial conditions which are arbitrarily close to an 
equilibrium where p pairs of neighbors are diametrically 
opposed, but for which convergence to a configuration where 
at most p — 1 pairs of neighbors are diametrically opposed 
is guaranteed; this reasoning implies that all configurations 
where neighbors are diametrically opposed are unstable, and 
therefore, only synchronized configurations are stable. A 
similar result on synchronization in SO( 3) is found in [20], 
and of synchronization on the sphere is found in [11], Also, 
Theorems 16 and [IT] do not provide any insight on whether 
the limits lim,^^ n ft) (for all i £ Af) exist; i.e, even if the 
unit vectors synchronize, it is not known whether those unit 
vectors converge to a constant unit vector or a time-varying 


unit vector. Some preliminaries results are found in XIV 


VI. Non-Tree Graphs 


In this section, we study the equilibria configuration in¬ 
duced by some more general, yet specific, network graphs. 
Also, we study the local stability properties of the synchro¬ 
nized configuration for arbitrary graphs. We first give the 
following definition. 

Definition 4: Given n vectors x ; G R 3 , for i £ {1, • • • , n}, 
we say that „ } belong to a common plane if there 

exists a unit vector u £ S 2 such that II (iz) x ; = X; for 
all i £ {l,--- , n}. We say that {x ; } ig{1 ... „ } belong to a 
common unique plane if there exists a single pair of unit 
vectors (+zz, —zz), with v £ S 2 such that II (zz) x ; = x £ for 
all i £ {1, • • • , n} 

Let us first discuss a property that is exploited later in this 
section. 


Example 3: Consider the distance functions d(ni,n 2 ) = 
/(l — rffn 2 ) where f(s) = a (|) for positive a and a > 1. 
For these functions, d max = a, f £ V°, and / G V 0 if a > 1 
and f £ Vo if a = 1. Suppose f k (s) = f(s) for all k £ M, 
and for some a and a. It follows that if r := < p 

(for some p £ M) and 


n(0) £ C 



arccos 



then Theorem’s flTIconditions are fulfilled and its conclusions 
follow. Similarly to example [T] by increasing a convergence 
for arbitrary initial values of rotational kinetic energy can 
be guaranteed; on the other hand, by increasing a we can 
increase the size of the cone where the agents need to initially 
be contained (up to a cone C (§))■ The main difference 
between this function and that in Example [T] is that, in this 
example, f £ V°. 


Under Theorem’s [FT] conditions, the group of agents can 
converge to configurations where one or more pairs of 
neighbors are diametrically opposed. However, it does not 
provide any insight on whether these equilibrium configura¬ 
tions are stable or unstable. For tree graphs, a configuration 
where p pairs of neighbors are diametrically opposed can be 


Proposition 18: Consider rii,n 2 G S 2 . If N(n 1 )n 2 7^ 0, 
then rn and n 2 belong a common unique plane. 

Proof: Consider iz = ll ^ ni | n3 ll £ S 2 , which is well 
defined since S (rq) n 2 7^ 0. It follows that n (iz) = in 
and that n (iz) n 2 = n 2 , which implies that rq and n 2 belong 
a common plane. Moreover, in and n 2 belong a common 
unique plane, since in and n 2 span a two dimensional space. 


Proposition 19: Consider in, ..., n„ G S 2 , with 
n' r n i+1 1 7^ 1 for all i = {1, • • • , n — 1}. If 

S (in) n 2 _ S (in,.,) n„ 

l!S(m)n 2 |r 

then all unit vectors belong to a common unique plane. 


Proof: 
l<n 3 | / 

\\S (n 2 )n 3 || 


Consider n = 3. Since |nfn 2 | 7^ 
1, it follows that || < S(n 1 )n 2 || 7^ 

^ 0. Additionally, by assumption, 

S (in) n 2 S (n 2 ) n 3 


±-, 


S (in) n 2 


= ±- 


|<5 (n 2 ) n 3 


1 and 
0 and 


(VI. 1) 


is satisfied. Consider then iz = ‘V ni ) n 2 

||o(ni)n 2 || 


G S 2 . It follows 


immediately that n (iz) in = n and that n(iz)n 2 = n 2 . 
It also follows that n (zz) n 3 = (I — iziz T ) n 3 = n 3 — 
















i'(v T n 3 ) = n 3 , where u T n 3 = 0 follows from taking 
the inner product of ( |VI| ) with n 3 . Altogether, it follows 
that n, n 2 and n 3 belong to a common unique plane (see 
Proposition [T~8j). For n > 3, it suffices to apply the previous 
argument n — 2 times. ■ 


Proposition 20: Consider n,, n„ £ S 2 . If 

±e 1 (n 1 ,n 3 ) = ••• = ±e„_ 1 (n„_ 1 , n„) then all unit 
vectors belong to a common plane, which is unique if 

ie 1 (n 1 , n 2 ) = ■ • • = ±e„_ 1 (n„_ 1 , n„) f 0 . 

Proof: If ±e 1 (n 1 , n 2 ) = ••■ = ±e„_ 1 (n„_ 1 , n„) f 0, 


it suffices to invoke Proposition 19 If ±e 1 (n 1 , n 2 ) = • • • = 
±e„_ 1 (n„_ 1 , n„) = 0 , it follows that irq = ••• = ±n„, 
and therefore, all unit vectors belong to a common plane. ■ 

Theorem 21: Consider the vector field ([HI]), the control 
law ( ]IV-E| ) with f k £ V° for all k £ M, and a trajectory 
x(-) of x(f) = f x (f, x(t), T o! (f, x(f))). If the network graph 
contains only independent cycles, then for each cycle, all its 
unit vectors converge to a common plane. 


Proof: By invoking Proposition 14 it follows that 
lim t _ > 00 (B®I)e(n(f)) = 0, i.e., e(n(-)) converges to the null 
space of B ®I. Now, consider a graph with only independent 
cycles and recall Proposition [4] Without loss of generality, 
consider that there is only one independent cycle and that 
the first n > 3 edges form that cycle. From Proposition [4] 
it follows that e(n) £ A f{B <g> I) =>■ ±e 1 ( 1 n, jii) = • •■ = 
±e„(„n, B n). In turn, from Proposition [20| it follows that all 
unit vectors that form the cycle belong to a common plane 
when (B ® I)e(-) =0. ■ 


Figure [7(a)1 exemplifies the statement in Theorem 21 where 
three agents form a complete graph, and thus there is one 
independent cycle, and where the distance function is the 
same for all edges. In this scenario, and because the distance 
function is the same for all edges, the unit vectors form an 
equilateral triangle. In Fig. [7Tb) four agents form a complete 
graph, which does not fit the conditions of Theorem 21 (a 
complete graph with fours agents has three cycles, but they 
all share edges with each other, i.e., they are not indepen¬ 
dent); for this scenario, we can find equilibria configurations 
where the unit vectors do not belong to a common plane, and, 


in Fig. 7(b) an equilibrium configuration is shown where the 
four agents form a tetrahedron. 




Fig. 7: Equilibrium configurations in a complete graph (1C 3 
in Fig. 7(a) and 1C 4 in Fig. |7(b)[ >, with the same distance 
function for all edges (in Example [3j with a = 1) 


Theorem 22: Consider the vector field ( |lll| ), the control 
law ( ]IV-E| ) with f k £ 'P" for all k £ A4, and a trajectory 
x(-) of x(t) = f x (t,x(t),T c ‘(t,x(t))). If the network graph 
contains only independent cycles or cycles that share only 
one edge, then all unit vectors belonging to each independent 
cycle converge to a common plane, and all unit vectors 
belonging to each pair of cycles that share only one edge 
also converge to a common plane. 


Proof: For network graphs with only independent 
cycles, the proof follows the same steps as those in the proof 
of Theorem [3T] Thus, and w.l.o.g, consider graphs with only 
two cycles that share only one edge, and recall Proposition [5] 
(adding other cycles does not affect the following conclu¬ 
sions, as long as those cycles are formed by other edges other 
than those that form the two previously mentioned cycles). 
W.l.o.g, assume that the first n = q+s— 1 > 5 edges are part 
of the two cycles that share only one edge, where the edges 
{!)■■■ 7 < 7 — 1, q} form the first cycle (|{1, • • • ,q-l,q}\ = q), 
while the edges {< 7 , q + 1 , • • • , n} form the second cycle 
f |{< 7 , g + l, ■ ■ • , 71 }| = s), and with the q th edge as the shared 
edge between cycles (for brevity, denote p = q — 1 and 
r = 9 + 1). Under the Theorem’s conditions, and by invoking 
Proposition 14 it follows that lim i _ >oc (.B ® I)e(n(f)) = 0, 
i.e., e(n(-)) converges to the null space of B(& I. From CxTI) , 
A f{B ® I) is spanned by the vector 

[(l,-i ® a) T (a + (3) T (l s _i ® f3) T O 1 ’] 3 ’, 

for any a,/3 £ R 3 (see Proposition [5]). 


Suppose „ 11 f ± ? n, where , ( n and ,,n are the unit vectors 
that form edge < 7 ; then e, ( („n, ,-n) f 0 and therefore a+(3 f 
0=>a^0V/3/0; moreover, q n and ,.n form a common 
unique plane (see Proposition [l 8 |. 


If a = 0 (or (3 = 0), it follows that ie^n, in) = 
••• = ±e p ( p n, p n) = 0 , and therefore the unit vectors 
in the first (or second) cycle belong to a common plane 
(Proposition 20 1 . Also, if a. = 0 (or (3 = 0), it follows 
that (3 f 0 (or a f 0 ), and therefore ±e r ( r n, f n) = 
••• = ±e„(„n, s n) f 0 , and therefore all unit vectors 
in the second (or first) cycle belong to a common unique 
plane (Proposition |20[ >. However, since the unit vectors of 
the first (or second) cycle are either all synchronized or some 
are diametrically opposed to others, it follows that all unit 
vectors in both cycles belong to a common unique plane. 


Consider now the case where a. f 0 and fl f 0, which 
implies that ie^n, in) = ••• = ±e p ( p n, p n) f 0 as 
well as ±e P ( P n, P n) = ••• = ±e„(„n, ft n) f 0. Thus, by 
Proposition [20] it follows that all unit vectors in each cycle 
belong to a common plane. Since , ( n and , ( n form a common 
unique plane, and since , ( n and , ( n belong simultaneously to 
both cycles, it follows that all unit vectors in both cycles 
belong to a common unique plane. 

Suppose now that 9 n = ± 5 n. Then e 9 ( 5 n, 9 n) = 0, and it 
follows that a + /3 = 0^>a = —/3 / 0Vq = (3 = 0. 
If a = (3 = 0, then e(n) = 0 <^> irq = • • • = ±n N and 

























the Theorem’s conclusions follow. If a = —{3 f 0, then 
±e 1 ( 1 n, in) = ••• = ±e p ( p n, p n) = ±e r ( r n, P n) = ••• = 
±e„(„n, „n) / 0, which implies that all unit vectors in both 
cycles belong to a common unique plane (Proposition [20|. 


Figure [8] exemplifies the statement in Theorem [22] with a 
network of six agents, with the network graph in Fig. |8(c)| 
and where the distance function is the same for all edges 
(distance function in Example [3] with a = 1). In this scenario, 
there are two cycles that share only one edge, one cycle 
composed of unit vectors {n-,, n 2 , n 3 , n 4 , n 5 }, a second cycle 
composed of unit vectors {n x , n 5 , n 6 }, and where the shared 
edge is formed by {rq, n 5 }. There are at least two equilibria 
configurations (apart from configurations where n, = ± n, 
for some i and j), which are given in Fig.|8(a)|and Fig. 8(b) 


where in both cases all unit vectors belong to a common 
plane. 


Propositions 21 and 22 focus on equilibria for some gen¬ 
eral, yet specific, network graphs. For arbitrary graphs, we 
can find equilibria configurations as exemplified in Fig. |7(b)| 
where four agents in a complete graph are at equilibrium 
when forming a tetrahedron. 


We now present a proposition, which will be useful in 
guaranteeing local asymptotic stability of incomplete attitude 
synchronization for arbitrary graphs. 

Proposition 23: Consider n £ C(a), for some a £ [0, |), 
and assume that i) the network graph is connected; ii) e(n) £ 
Af{B 0 I), with e(-) as defined in ( |IV-B| >. This takes place 
if and only if 3v £ S 2 : n = (1 N <g> v). 

Proof: For the sufficiency statement, it follows that, if 
3u £ S 2 : n = (1 N 0 v), then all unit vectors are contained 
in a |-cone, i.e., n £ C(|); and, moreover, e(l w 0 is) = 
0 C A f(B 0 I). 

For the necessity statement, the proof is as follows. For a 
tree graph, (B 0 I)e(n) = 0 e(n) = 0 follows. This 
implies that 10 = in, for all (i,j) £ Af x A f, but since 
n £ C(|), it follows that n, = rr, for all (i,j) £ Af x Af t . In 
a connected graph, this implies that n, = n, for all (i,j) £ 
Af x Af, and therefore 3u £ S 2 : n = (1 N 0 u). 

For an arbitrary graph, the null space of (B 0 I) may 
be more than 0, i.e., (B 0 I)e(-) = 0 e(-) = 0. We 

anticipate the final result by stating that if n £ C(^), then 
(.B 0 I)e(n) = 0 => e(n) = 0, in which case we conclude 
again that 3u £ S 2 : n = (l w 0 u). Consider then an 
n £ ( S 2 ) N , such that (B 0l)e(n) = 0. This means that, for 
every i £ Af (B i: stands for the i th row of B), 

0 = (B i: ®I)e(n) ^S(n t ) £ ^ nfn^nj.(VI.2) 

Since n £ C(a), it follows that there exists a unit vector 
H £ S 2 , such that /i r n ; > cos(a) > 0 for all i £ Af. 
Taking the inner product of ( |VI[ i with S (n,) /r„ it follows 
that /x T II(n,)J] jG jv i /E,j)( 1_1 E n j)n,- = 0, which can be 


expanded into 

- n X) (mX - (M r n 1 )nfn J ) = 0(.VI.3) 

Now, consider the set T = {i £ Af : i = argmax ie j\/(l — 
/T'n,)}, and choose k £ T (in the end, we show that, in 
fact, T = Af). Notice that 0 < cos(a) < p T n k < /i T n ; for 
all k £ T and all j £ Af. As such, it follows from ( |VI| ) with 
i = k that 


0 < cos(a) J2 - nfn,) (1 -<n,) < 

jeU k 

< E fUi.A 1 ~ n fw) (M r n J - (/i T n fc Xn j ) = 0(VI.4) 

oeM k 

Notice that the lower bound (on the left side of m is 
non-negative and zero if and only if all neighbors of agent 
k are synchronized with itself (note that lim s _ >2 - 
may be 0, but since n £ C(a), /' (i i) (s)| s _i_ n T n . can only 
vanish if s —> 0 + ). As such, it follows from ( |VI[ > that all 
neighbors of agent k are contained in T, i.e., Af k C T. As 
such, the same procedure as before can be followed for the 
neighbors of agent k, to conclude that the neighbors of the 
neighbors of agent k are all synchronized. In a connected 
graph, by applying the previous reasoning at most N — 1 
times, it follows that all unit vectors are synchronized, or, 
equivalently, that 3iz £ S 2 : n = (1 N 0 u). ■ 


Proposition [23] has the following interpretation. Recall that 
{x £ ( S 2 ) N x R 3JV : ( B 0 I)e(n) = Cfu;, = 0 for i £ 
£.11(10) w, = 0 for j £ £} C 12^, where f2^ is the set 
of equilibrium points. For example, we have seen that, for 
specific graphs, all equilibrium configurations are such that 


all unit vectors belong to a common plane (see Theorems 21 
and|22|, as illustrated in Fig. [8] However, if we can guarantee 
that along a trajectory x(-) of x(t) = f x (t, x(t), T cl (t, x(t))), 
3a £ [0, |) : n(t) £ C(a),\/t > 0, i.e., if we can guarantee 
that all unit vectors remain in an closed a-cone smaller 
than an open f-cone, then we can invoke Proposition 23 to 
conclude that lim t _ >oo (i30l)e(n(f)) = 0 => liinj^^jnAf) — 
10(f)) = 0; i.e., that convergence of e(n(-)) to the null space 
of 5 01 implies synchronization of the agents. 


This motivates us to introduce a distance d* > 0, which 
will be useful in guaranteeing that, along a trajectory x(-), 
3a £ [0, |) : 11(f) £ C(a),\/t > 0. Consider then 

d * = min/„ (l - cos < d min , (VI.5) 

which satisfies ,/jr 1 (d‘) < 1 — cos ^ aEl) f° r ^ £ AT 
Notice that d* < d min , since d min = min fe£jV , lim a _ >2 - f k {s), 
1 — cos( 3 ) < 2 for all N > 2, and since all f k (-) are 

increasing functions in (0,2). As shown next, if .D(n(0)) < 
jj (and H(uj(0) = 0), then the network of unit vectors is 
forever contained in a closed a-cone, for some a £ [0, ((). 

Theorem 24: Consider an arbitrary network graph, the 
vector field the control law ( |IV-F[ i, and a trajectory 

x(-) of x(t) = f x (f, x(f), T cl (f, x(f))) If x(0) £0.1 = 
{x £ O^ x W N : V(x) < d*} then synchronization is 
asymptotically reached, i.e., lim t _ >ao (n i (f) — n, (f)) = 0, for 
all (i, j) £ Af 2 . Moreover, all implications of Proposition [To| 
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(a) Equilibrium configuration with 
network graph shown in Fig. |8(c)| 


(b) Equilibrium configuration with 
network graph shown in Fig. |8(c)| 


(c) Graph with two cycles that 
share only one edge 


Fig. 8: Two equilibrium configurations for group with network graph shown in Fig. 8(c) with the same distance function 

for all edges (in Example [3] with a = 1) 


also follow. 


VII. Simulations 


Proof: Since d * < we can invoke Proposition [To] 
and infer that lim t _ >00 (f? (g) I)e(f) = 0 (as well as all other 
implications stated in the Proposition). Since V"(x(-)) < 0, 
it follows that f k ( 1 - fc n T (-) s n(-)) < D{ n(-)) < V(x(-)) < 
1/(0) < d*, for all k £ Ai. In turn, this implies that 
6>( fc n(-), s n(-)) < arccos(l - /* -1 (d*)) < fjvti- for all 
k £ Ai. Since the angular displacement between any 
two unit vectors n and n, in a connected graph satis¬ 
fies ffjnf-) ,n.j(-)) < (IV- l)max feeM 0( fe n(-)> *“(•)) (see 
Proposition 27i, it follows that sup t>0 6>(n, (f), (£)) < | 


between any two unit vectors n, and n . As such, it follows 
from Proposition ul that n(-) £ C(| sup t>0 rr,(f))), 

where } sup (>0 f/(ri,(f), n,(f)) < Finally, we invoke 
Proposition [23] which implies that (B Q I) e (n(f)) = 

0 => lim t _ foo (n i (f) — !!,■(£)) = 0. ■ 


Let us provide a corollary to Theorem 24 with an easier 
to visualize region of attraction. 


Corollary 25: Theorem 


if 


24 


holds if r := < 1 and 


n (0) 6 C arccos ^1 - min f k 1 (1 - r )))) 

with d* as in m 

For proving Corollary [25] it suffices to check that if its 
conditions are satisfied, then V r (x(0)) < d*. 

Remark 26: Comparing Theorems [l6| an d [24] it follows 
that the region of attraction in Theorem |l6| is larger than that 
in Theorem |24| Loosely speaking, the region of attraction in 


Theorem 


16 


attraction in Theorem 
network graph topology, and in fact, a tree network graph 
provides stronger results. 


Theorems [16] and 24 provide asymptotic results, namely 
lirn,^ e(n(f)) = 0. Remark 34 in Appendix provides some 


We now present simulations that illustrate some of the 
results presented previously. For the simulations, we have a 
group of ten agents, whose network graph is that presented 
in Fig. im The moments of inertia were generated by adding 
a random symmetric matrix (between —I and I, entry-wise) 
to the identity matrix. For the initial conditions, we have 
chosen u>(0) = 0 and we have randomly generated one 
set of 10 rotation matrices. In Fig. [9] h, = [10 0] T for 
i = {1,2,3,4,5} and h, = [010] T for i = {6,7,8,9,10}, 
and since these are not necessarily principal axes, we apply 


the control law ( |IV-L[ i, with £ = 0 and £ = Af. In Fig. 10 


rq is the principal axis of , with largest eigenvalue, for 
i = {1,2,3,4,5}, and n, = [10 0] r for i = {6,7,8,9,10}. 
Therefore, we apply the control law ( |IV-E[ ), with £ = 
{1,2,3,4,5} and £ = {6,7,8,9,10}. For the edge 1, 
we have chosen f 1 (s) = 10 tan 2 (0.5 arccos(l — s)). For 
the other edges, we have chosen f k (s) = 5s, for k = 
A4\{1}. Notice that we have chosen a distance function (for 
edge 1) that grows unbounded when two unit vectors are 
diametrically opposed. As such, it follows that agents 1 and 
6 will never be diametrically opposed, under the condition 
that they are not initially diametrically opposed. We have 


also chosen <r(x) = k- 


with k = 10 and <j x = 1. 


is ^p- > 1 times larger than the region of 
This difference comes from the 


insight on exponential convergence to 0. 


y/ 

For this choice, we find that er max = ka x = 10. As such, 
for all agents, except 1 and 6, an upper bound on the norm 
of their torque is given by a ln;ix + 2 • 5 = 20 (the factor 
2 relates to the fact that all agents, except 1 and 6, have 
two neighbors, and the factor 5 comes from f k (s) = 5 s => 
f' k (s ) = 5). For agents 1 and 6, no upper bound can be found 
(more precisely, a bound can be found, but it depends on 
the initial conditions). Given these choices, it follows from 
Corollary 25 that if n(0) £ C(~ 1°) then synchronization 
is guaranteed. We emphasize, nonetheless, that Corollary [25] 
provides conservative conditions for synchronization to be 
achieved, and the domain of attraction is in fact larger. We 
also emphasize that, for tree graphs, the domain of attraction 
is considerably larger: for example, if we removed the edges 
between agents 1 and 2, and between agents 6 and 7, we 
would obtain a tree graph, and Corollary [13] would read 







































(c) Error angle between 
neighbors without noise 


(d) Error angle between 
neighbors with noise 


(c) Error angle between 
neighbors without noise 


(d) Error angle between 
neighbors with noise 


Fig. 9: Synchronization in network of 10 unit vectors with 
and without noise, where blue agents perform 
synchronization of their first axes (h, = [10 0] T ) and black 
agents synchronization of their second axes (h, = [010] T ). 


as n(0) £ C(~ 18°). Finally, we emphasize that we can 
increase the size of the cones in Corollaries [13] and [25] 
by choosing different distance functions, as exemplified in 
Example |T] 


For each Fig. [9] and Fig. 10 we provide two simulations: 
one where the control law is that in ( |IV-E| i and another where 
the control law is that in ( |IV-E| ) but corrupted by noise; 
specifically, for each agent i £ J\f, T 4 (f) = Tf(t,x(t)) + 
0.1AJ0 0 1] T , where A ; corresponds to the largest eigenvalue 
of J,. 


The trajectories of the unit vectors for described conditions 
are presented in Figs. 9(a)-9(b) and 10(a) - 10(b) Notice that 


despite not satisfying conditions of Theorem 24 (the unit 
vectors are not always in a | cone), incomplete attitude 
synchronization is still achieved. This can be verified in 
Figs. |9(c)H9(d)| and |10(c)|[T0(d)] which present the angular 
error between neighboring agents. In Figs. 10(a) and |10(b~ 


the control laws are different between agents 1-5 and 6-10. 
The former perform synchronization of principal axes, by 
applying the constrained control law ( |IV-E| i; while the later 
perform synchronization of their first axes, i.e., n, = [10 0] T , 
by applying the control law ( |IV-C| ). In Figs. 9(d) and 1 10(d)] 
for which the control laws were corrupted by noise, perfect 
synchronization is not asymptotically achieved. Instead, the 
unit vectors converge to a configuration where they remain 
close to each other (error of « 5° between neighbors). 


VIII. Conclusions 

In this paper, we proposed a distributed control strategy 
that guarantees attitude synchronization of unit vectors, rep¬ 
resenting a specific body direction of a rigid body. The 
proposed control torque laws depend on distance functions 


Fig. 10: Synchronization in network of 10 unit vectors with 
and without noise, where blue agents perform 
synchronization of principal axes and black agents 
synchronization of their first axes (i.e., h = [10 0] T ). 



Fig. 11: Graph with 10 agents, where edge 1 is formed by 
agents 1 and 6 

in S 2 , and we provide conditions on these distance functions 
that guarantee that i) a synchronized network is locally 
asymptotically stable in an arbitrary connected undirected 
graph network; ii) a synchronized network can be achieved 
for almost all initial conditions in a tree graph network. We 
imposed conditions on the distance functions that guarantee 
that these are invariant to rotation of their arguments, which 
means that the proposed control laws can be implemented by 
each individual rigid body in the absence of a global common 
orientation frame, i.e., by using only local information. Addi¬ 
tionally, if the direction to be synchronized is a principal axis 
of the rigid body, we proposed a control law that does not 
require full torque actuation, and, more specifically, it only 
requires torque in the plane orthogonal to the principal axis. 
We also studied the equilibria configurations that come with 
certain types of graph networks. Directions for future work 
include studying the stability of all equilibria configurations, 
apart from the synchronized configuration; to determine 
whether a synchronized network converges to a constant unit 
vector in a fixed, though unknown, orientation frame; and to 
extend the results to complete attitude synchronization. For 
complete attitude synchronization, though, it is not possible 
to construct a constrained control law, since it requires full 
torque actuation. 
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IX. Triangular Inequality 


Proposition 27: Consider three unit vectors n , n 2 and n. 
The following triangular inequality is satisfied, 

0(ni, n 2 ) < ^(n!, n) + 0(n, n 2 ), 

where 0 : S 2 x S 2 H > [0,7r], defined as 0(ni,n 2 ) = 
arccos(nfn 2 ). 


Proof: Since n and n 2 can be written as 
lii = cos(0(n!, n))n + sin(0(n!, n))u u 
n 2 = cos(0(n 2 , n))n + sin(0(n 2 , n ))u 2 , 

where v ± , u 2 £ S 2 are unit vectors orthogonal to n, then 


= 008 ( 0 ( 11 !, n)) cos(0(n 2 , n)) + sin(0(n!, n)) sin(0(n 2 . n))i/^i/ 2 
= cos(0(n l5 n) + 0(n 2 , n)) + sin(0(n 1 , n)) sin(0(n 2 , n)) (1 + vfi/f) 


> cos(0(n!, n) + 0(n 2 , n)). 

As such, ^(n!, n 2 ) < 0(n 1; n) + 0(n 2 , n). ■ 

Proof: [of Proposition |7) Let us construct an u € S 2 
such that 0( n k ,i/) < a for all k € Af , with 0(-,-) defined 
in Proposition 27 If max (j j)6Af 2 0(ru, n 2 ) = 0, it suffices 
to pick v = n , for any i £ Af, and the Proposition’s con¬ 
clusion follows. If 0 < max (i j)EAf ! 0(n i ,n J ), take (k,l) = 
argmax (i 0(n„n J ), and consider then u = 
which is well defined since 0 < 0(n fc ,n i ) < |7T < 7r. By 
construction, it follows that 9(n k ,u) = 0( n t ,v>) = 

From Proposition [27] it follows that, for any j £ Af, 
0( n !, v) < 0(n j ,n fc )+0(n fc ,i/) < 0(11;, n fc ) + 0(n fc , u) < 
§0(n fe ,n,). Thus, if 0(n fc ,n,) < | a, for some a £ [0,7r], 
then 0(11.. u) < a for all k £ Af, and the Proposition’s 
conclusion follows. ■ 


X. Solution to PDE (IV. 1) 

The partial differential equation ( |IV-B| ) is a restriction 
imposed on the distance functions. The question that follows 
is what kinds of distance functions d : S 2 x S 2 —» R+ 
satisfy (|IV-B|). Distance functions of the type 1/(11,. n 2 ) = 










/(I - nfn 2 ) (with / G C 1 ((0,2),K> o )) satisfy ( |IV-B| ). 
In fact, those are the only type of distance functions that 
satisfy ( ]IV-B| ), as we show next. 

Proposition 28: Consider two orthogonal unit vectors, i.e., 
(v 1 ; v 2 ) G n = {(rii,n 2 ) G S 2 x S 2 : n}n 2 = 0}. 
Additionally, consider a function 77 G C 3 (R 3 x R 3 ,R >0 ) that 
satisfies 

S ( v,)*^ + S(v,)*?!^ = 0, (X.1, 

ov 1 av 2 

for all (v 1 ? v 2 ) E Then, rj(v u v 2 ) is constant for all 

(vd v 2 ) g n. 

Proof: The function ?y(-, •) is constant for all (v,, v 2 ) G 
Q, if the differential dpf,-) is zero along (v!,v 2 ) G Cl. 
This condition is proved next. Taking the inner product 
of {28} with V; (i G {1,2} and j = {1,2}\{*}), it follows 
that vfS (v,) ^’ V2) = 0, which means dvi g^ V2) is 
orthogonal to both iS(v 1 )v 2 and ( S(v 2 )v 1 . Since v l5 v 2 
and S (vj v 2 form a basis of R 3 , it then follows that 

dri(v i,v 2 ) _ ( T dr 1 (v 1 ,v 2 ) \ , ( T 8 i)(vi,v 2 ) \ 

avi — } v i a Vl { Vi “r } v 2 a Vl ) 

=/i 1 (v 1 ,v 2 )v 1 +/i 3 (v 1 ,v 2 )v 2 , (X.2) 

a?;(vi,v2) _ ( t aT)(vi,v 2 ) \ , ( T a»?(vi,v 2 ) \ 

av 2 ~ } v i av 2 { 1 v 2 av 2 ) v 2 

= /i 4 (v 1 ,v 2 )v 1 + /t 2 (v 1 ,v 2 )v 2 . (X.3) 

By replacing © and © in ( [28] ), it follows that 
h 3 {v i,v 2 ) = (v 2 , v 2 ) = /i(v 1 ,v 2 ) for all (v 1; v 2 ) G 

(notice that <S( v i ) v 2 0). Now, recall that (v l7 v 2 ) G fi 
means that vfvi = 1, vjv 2 = 1 and v}v 2 = 0. 

In turn, this implies that dvfv 1 = 0, dv 2 v 2 = 0 

and dv}v 2 + dv 2 v 1 = 0. From ©, it follows 
that dvf 9t?( 3 ^ V2) = h 1 {-,-)dvfst 1 + /i(-,-)dv}v 2 = 

h(-, -)dvfv 2 , and, from ©, it follows that rfvj = 

h(-,-)dv 2 v 1 + h 2 f, -)rfv}v 2 = hf,-)d'v 1 2 -v 1 . Since 

d*?( Vl ,v a ) = dv- d -^l + 

/i(vi, v 2 ) (dv[v 2 + dvjvi) = 0, the Proposition’s 

conclusion follows. ■ 


Proposition 29: Consider a function 77 G C X (R 3 x R 3 , R >0 ) 
that satisfies 

c* \ ^ 7 ?( n 1 , n 2 ) c \ ^(ni, n 2 ) „ , v 

S ( n i)- q-1- s ( n 2 -x-= 0, (X.4) 

OTj a n 2 

for all (n!,n 2 ) G S 2 x S 2 . Then, 77 (n 1 ; n 2 ) = h( 1 — n}n 2 ) 
for (n l5 n 2 ) G S 2 x S 2 and for any h G C 1 ((0, 2),R> 0 ). 


Proof: Condition ([29} is satisfied for (n,, n,) G S ! x 
<S 2 , and, in particular, it is satisfied for (n!,n 2 ) G 0 = 
{(vi, v 2 ) G S 2 x S 2 : (v}v 2 ) 2 / 1} which is a set of 
full measure w.r.t. S 2 x <S 2 . For (n,, n 2 ) G 0, consider the 
variables 

n, + n 2 n 2 - n 2 T 

v, = x---x, v 2 = tt -= n, n 2 , 


Hrii + n 2 1 | Uni — n 2 || 

and notice that, by construction, ru and n 2 can be obtained 
from v 1; v 2 and v 3 (indeed, n 12 = 83 V! ± \J\ — v 3 v 2 ), 
meaning that a bijection exists between (rii,n 2 ) and 
(v 1; v 2 , v 3 ); moreover, v 2 is orthogonal to v 2 , i.e., vfv 2 = 0 . 


The proof hereafter explores the previous change of vari¬ 
ables. First, notice that there exists a function h(-, •, •) such 
that 77 ( 11 !, n 2 ) = /i(v 1 (n 1 ,n 2 ),v 2 (n 1 ,n 2 ),?; 3 (ni,n 2 )) or al¬ 
ternatively 77 ( 11 !(v l7 v 2 , v 3 ), n 2 (nj, n 2 , v 3 )) = h(vi, v 2 , v 3 ). 
Then, from ([29}, it follows that (for simplicity, all arguments 
in the equations below are omitted) 


5 (nf) 


S (n 2 



dv 1 dh dv 2 dh 
d n}<9v{ 

3v, dh 
dn 2 dv 1 


= 0 (X.5) 


e l sw st tsw ik 


S (n,) —- + S (n 2 ) —- —— 


dh 

{S (n,) n 2 + S (n 2 ) nj — = 0 

_, . dh _, . dh 

&S (v x ) --h S (v 2 ) — = 0 

av! av. 


(X. 6 ) 


(X.7) 


where the chain rule has been applied to derive © 
from {29}; and from {X} to ©. the relations (recall that 

& (nirj = ( : _ mm) m } 

dv, , 1 , 

— = (I — VjV ) TT- ——tt, fori G {1,2}, 

d n, ||ri! -f n 2 || 

X- 1 = ± (I - v 2 v 2 t ) r- --- TT, for j G {1,2}, 

on, || n, — n 2 1| 

together with the relations S (ii,) v 2 + S ( n,) v = 0 and 
S (ii,) v 2 — S (n 2 ) v 2 = 0 have been used. It follows from 
Proposition 28 and from ( fx) . that /v.( v,. v 2 , v 3 ) does not de¬ 
pend on v x and v 2 , and, as such, 77 ( 11 !, n 2 ) = h(-, •, n}n 2 ) = 
h(l — nfn 2 ). ■ 


XI. Auxiliary Examples (Section |IV-A} 


Example 4: In Fig. 12 three graphs with equivalent in¬ 
cidence matrices are presented (incidence matrices in (|4}, 
0 and {f}). (Two matrices D,E G R ^ 171 are said to 
be equivalent iff there exist non-singular P G R" x " and 
Q G R m * m such that D = PEQ [27].) The incidence matrix 
in Fig. 12(a) is given by 


B 1 


1 0 1 
-110 
0 -1 -1 


(XI. 1) 


Consider a graph, with incidence matrix B 1 G R NxM , and a 
second graph, with incidence matrix B 2 G R NxM , obtained 
from the previous after a swap of nodes. Then If = T, /i 2 I, 
where F, G R JVXJV is a permutation matrix corresponding to 
the swap of nodes (permutation of rows) and I G R MxM is a 
diagonal matrix with diagonal elements ±1 (column k G A4 
needs to be multiplied by —1 when, after permutation of 
rows, i becomes larger than j with (i,j) = tc - 1 (fc)). The 
graph in Fig. |12(b)| is obtained from that in Fig. 1 12(a) | by 
swapping nodes 1 and 2 , and, as such, its incidence matrix 







































is given by 


B 2 =1^1 = 


0 1 0 
1 0 0 
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-1 


(XI.2) 

where column 1 needs to be multiplied by -1 since rows 1 
2 have permuted and k( 1,2) = 1. Consider again a 


and 


graph, with incidence matrix B, £ R N x M , and a second 
graph, with incidence matrix B 2 £ R NxM , obtained from 
the previous after a swap of edges (edges’ numbers). Then 
B 1 = B 2 T 2 , where T 2 £ R M x M is a permutation matrix cor¬ 
responding to the swap of edges (permutation of columns). 


The graph in Fig. 12(c) is obtained from that in Fig. 12(b) by 
swapping edges 1 and 3, and, as such, its incidence matrix 
is given by 


B , = -BoT, = 


1 

1 

o' 


0 

0 

T 


0 

1 

T 

-1 

0 

1 


0 

1 

0 

= 

1 

0 

-K 

0 

-1 

-1 


1 

0 

0 


-1 

-1 

0 


Example 5: Figure [3(a)| displays a graph with two indepen¬ 
dent cycles, with incidence matrix and its null space in 0- 


B 1 = 


10 1110 
-1 1 0 0 0 0 

0 - 1-1 0 0 0 
0 0 0 -1 0 1 

0 0 0 0 -1 -1 


Af(B x ) = span 


xXI.4) 


Example 6: Figure |3(b)| displays a graph with one inde¬ 
pendent cycle and two cycles that share only one edge, with 


1 

o ••• 

0 

1 ■ 

-1 

t ■■■ 

0 

0 

0 

-l 

0 

0 

0 

o ••• 

1 

0 

0 

o ••• 

-1 

-1 

0 

0 0 

0 

0 


Be = 0. Without loss of generality, assume that 


B c ±[B' e j Bl ] = 


which means, IP = EEi ~ 1 (B 1 c ) i . Consequently, the non¬ 
zero vector e = [0^,_ M 1^ _ 1 — 1] T belongs to the null- 
space of B , since [B e Bl Bf]e = Ei'=i _1 (-®c)i — Bl = 0. 
For B 0 I„, its null-space is spanned by the columns of 
e0>I n = [0I n • • • I ra — l n ] T . In general, the incidence matrix 
will be of the form B = T, [B, B) Bf\\T 2 , with T, £ R NXJV 
and r 2 £ R MXM as permutation matrices and I £ ]R AfxM as 
a diagonal matrix with diagonal elements ±T, i.e., the inci¬ 
dence matrix B will be equivalent to [B r Vi) Bf\ apart from a 
reordering of nodes’ and edges’ numbering. In that case, the 
null space of B will be spanned by T^I[0^_ Mc l^ c _ 1 — 1] T , 
or, alternatively, Af(B) = {e £ M M : e k = ±e,,Vfc, ( £ C}, 
where C is the set of the edges that form the cycle. For m 
independent cycles, (C' i } i={1 .. m} , i.e., cycles that do not 
share edges, Af(B) = {e £ M M : e k = ±e t ,fk,l £ Ct,i = 
{1, • • • , m}}. Similar conclusions extend to Af(B 01). ■ 

Proof: [of Proposition [5j Without loss of generality, 
consider a graph with only one pair of cycles that share 
only one edge. Similarly to the proof of Proposition [4] we 
partition the incidence matrix as B = [B E B, and through 
a similar argument we conclude that |A/"(.B)| = 2. In this 
scenario, B, is decomposed in three parts. Without loss of 
generality, assume that 


B, = { B\ Bl j VP 


B 1 


1 

0 

0 

0 

0 

-1 


0 

0 

1 

0 

0 

-1 


Af(B) = span < 


10(5, 

defined in 0). 


■ 1 

-1 

0 ■■■ 

1 ••• 

0 

0 

0 

0 

1 

0 

0 ■■■ 

0 ••• 

0 

0 

1 ■ 
0 

1 


0 


0 












1 


0 


0 





-1 0 


0 

0 




-1 


1 


0 




0 

0 ■■■ 

1 

0 

0 

0 ••• 

0 

0 

0 


0 


1 

-1 


• ■ (XI.5) 


0 

0 0 

-1 

1 

0 

0 0 

0 

0 






— 

0 
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1 ••• 

0 
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-1 
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0 ■■■ 
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-1 ■■ 
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0 ■■■ 
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0 ■■■ 
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XII. Proofs of Propositions 


Proof: [of Proposition [4j Without loss of generality, 
consider a graph with only one cycle composed of M c edges. 
Then, its associated incidence matrix B can be divided in two 
parts, i.e., B = [B E B,]: B, corresponds to all edges that are 
part of the cycle; and B s corresponds to all other edges. B, 
can be further partitioned in two parts, i.e., B c = [Bl VP], 
where Bl is a single column matrix corresponding to a single 
edge of the cycle and such that [B s B) \ corresponds to a 
tree. Since [B E B) ] is a tree, it must be positive definite (see 
Proposition [2]); thus, it follows that the null space of B has 
dimension 1, i.e., \J\f(B) \ = 1, and there exists a unique (up 
to a scalar multiplication) non-zero vector e £ R M such that 


where {TP £? adj ] corresponds to one cycle (of dimension 
M c > 3), \1P B„il corresponds to the second cycle (of 
dimension M C2 > 3), and /f id| corresponds to the shared 
edge. As such, B adj = E ZT\ B Vp = which 

means that M{B) is given by 



(XII. 1) 


while A f(B 0 I rt ) is spanned by the columns of 


0 


0 

1m c , -i 0 In 


0m c1 -i ® In 

o 

s 

n 

o 

1 

1 

1m c „ — i & In 

-In 


-In 


(XII.2) 






























































(a) Graph with incidence (b) Same graph as in (c) Same graph as in 

matrix B\ in (XI. 1) Fig. 112(a)| apart from a swap Fig. 112(b)] apart from a swap 

between nodes 1 and 2 between edges 1 and 3 
(incidence matrix B 2 (incidence matrix B 3 

in (XI.2)). in (XI.3)). 


Fig. 12: Three graphs with equivalent incidence matrices. 


The null space in ( |X11| > can be equivalently written as 
[° T al u ci _i P 1 ^-! - (« + P)] T for an y € R. In 
general, the incidence matrix will be of the form B = 
I\ [B £ B\ B 2 o -Eyir 2 , with r! £ R NxJV and r 2 £ R MxM as 
permutation matrices and I £ R M x M as a diagonal matrix 
with diagonal elements ±1; i.e., the incidence matrix B 
will be equivalent to [71,-, B\ IP £? adj ] apart from a reordering 
of nodes’ and edges’ numbering. With that in mind, the 
null space in ( |XII| ) can also be written as Af(B) = {e £ 
R" : e fe = ±e,, Vfc, l £ C'\{C' n C 2 }, e p = ±e„Vp,q £ 
C' 2 \{C ,1 nC ,2 } ; } , with C 1 = {M-M C1 -M C2 + 1,..., M — 
M C2 - 1} U {M} and C 2 = {M - M C2 ,..., M - 1} U {M} 
(i.e, C 1 corresponds to the set of edges of one cycle, and C 2 
corresponds to the set of edges of the other cycle, where the 
M th edge is the shared edge). Including independent cycles 
or other pairs of cycles that share only one edge does not 
affect the previous conclusions regarding the null space, and, 
therefore, the conclusions of the Proposition follow. ■ 

XIII. Instability 

Under Theorem’s [TT] conditions, a configuration where two 
(or more) neighboring unit vectors are diametrically opposed, 
with all the other neighboring unit vectors synchronized, is 
an equilibrium solution. In fact, consider a constant unit 
vector n*, u> = 0 and n, = ±n* for all i £ Af. Under 
these circumstances, TJ = 0 for all * £ Af (see ( |IV-C| l), 
implying that all those configurations correspond to equi¬ 
librium configurations. Multiple equilibria exist because the 
control laws ( |I V-C[ > are continuous and the states n, (for all 
i £ Af) evolve in a non-contractible set (see [21]). Notice, 
however, that the equilibrium u> = 0 and n, = n* for all 
* £ AT is guaranteed to be (locally) asymptotically stable 
(see Theorem [24} . For all other equilibria, their stability, or 
lack thereof, is a subject of current research. 

In this section, we consider simpler dynamical agents and 
study the stability of equilibria where two (or more) neigh¬ 
boring unit vectors are diametrically opposed. Consider then 
agents n,, for i £ AT, with kinematics and whose angu¬ 
lar velocities ui are assumed to be controllable (intuitively, 
this corresponds to controlling single integrators rather than 
double integrators). Consider also the Lyapunov function 
U(n) = S (fc n , in), whose time derivative yields V = 
w T 7J(B®I)e £Li d k ( k n, in). This motivates the choice of 


the control law [u>f ... u)%] T = lj = —7 Z T (B Cg) I)e, since 
it implies that V = - || (B ® I)e|| 2 < -X laiI1 (B T B) ||e|| 2 . 
Since V{t) is uniformly continuous, Barbalat’s Lemma suf¬ 
fices to conclude that e converges to 0, which, in turn, 
implies that there exists a (possibly time-varying) unit vector 
n*(t) such that lim^oo (il (t) ± n*(f)) = 0, for all i £ Af. 

Let us focus on the equilibria configuration where n* is 
time invariant, i.e, focus on scenarios where all the unit 
vectors converge a constant unit vector, either +n* or —n* 
(the study of stability, or lack thereof, for time-varying iT (7) 
is a subject of current research). For the remaining part of 
this section, it is shown that the previous equilibrium config¬ 
uration is unstable whenever two (or more) neighboring unit 
vectors are diametrically opposed. This conclusion follows 
from a linearization procedure, which is described next. For 
a constant n*, consider the variables 

x, = n(n*)n,, (XIII. 1) 

Si = sign (nfn*), (XIII.2) 

for all i £ Af. From ( |XIII| ) and ( |X 11 [] >, it is possible to 

reconstruct rq, since 

n, = (nfn*)n* + II (n*) n, = (nfn*)n* + a^XIII.3) 
= sign (nfn*) ^/l — Ha^H 2 ]!* + x { (XIII.4) 

= Si\Jl- ll^fn* + Xi (XIII.5) 

where the relation (rLn,,) 2 + ||II (rii) n 2 1|“ = (njiia) 2 + 
||(S (nj) n 2 1| 2 = 1 has been used from ( |XIII| i to fjXTTTji. Since 
njn* can only change sign when n ; is orthogonal to n*, it 
follows that 

^ = 0, (XIII.6) 

if ||aJi|| < 1 for all time (||Xi|| = \J 1 — (n'/ n') 2 ). Addition¬ 
ally, from (JXlTTJi, it follows that 

nfn fc = ( Si yjl - \\Xi || 2 n* + a ',i) T {s k \J 1 - ||x t [| 2 n* + x k ) 



for all i,k £ Af, and where the fact that, for all i £ Af, x t is 
orthogonal to n* has been used. For x t = x k = 0, = 

sign ((nfn*)(n^n*)) = 1 if n, and are synchronized 
and SiS k = sign ((nf n*)(n^n*)) = —1 if n ; and n fc are 
diametrically opposed. Thus, by continuity, for small ||at { || 
and ||a; fc ||, SiS k = sign ((nfn*)(n^n*)) = 1 if n ; and 






























n fe are almost synchronized (i.e., 1 — n'n, is small ) and 
SiS k = sign ((nf n*)(n£n*)) = —1 if n, and n fc are almost 
diametrically opposed (i.e., —1 — n' n, is small). 


Since n* is constant, it follows from ( |XIII| ) that 

sc, =n (n*) fi, 

= n (n*) (-S (n,)w,) 

=n (n*) (-S (n 4 ) S (n<) EfceM n *)K) 

=11 (n*) | (I — n ; nf) ^ ff ik) (d e {n l7 n k ))n k 

= E /E,fc)( d »( n '’ n fc))( n ( n *) n k ^ ( n fn k )II (n*) n,) 

= ( x i - (nfn fc )a\) 

= ( x i - x <) + 

k^Mi 

fUi,k)( d »( n ^ n k)) (! - (nfn fc ))® 4 , (XIII.7) 

k^Mi 

where n and n fc can be expressed as in ( |XIII| ). By stack¬ 
ing all s, as a; = [tef ... as^.] and as s = 

[sf ... 4] T , (jXTTTji and ( |XIII| ) can be written compactly 

as 


x = — (B <8> 1)1?! (at, s)(B (g) I) T a; + D 2 (x 1 s)x 
=:A(x,s)x, (XIII.8) 

s =0 if H®,|| < 1 for all ieAT (XIII.9) 


with 

D 1 {x 1 s) 

D 2 (x,s) 


f' k (d e ( k n, s n)) (g)I 


fUiA d °( n *’ n i))0--( n T n k)) ) 

i€JVi / 


where both I), (x, s) and l) 2 (x. s) are positive semi-definite. 
As expected, the states x* = 0 and s* = 1 £ {z £ R w : 
z? = 1, V* = {1, • ■ ■ , N}} are equilibria of the system with 
dynamics ( |XIII| ) and ( |XIII| ). The linearized system for ( |XIII[ ) 
and (|XIII|), and around each equilibrium, is given by 


X 


A(0, s*) o' 

X 

s 


0 0 

S 


(B®I)Z? 1 (0,s*)(B(g>I) T 

0 


(XIII. 10) 


z? 2 (0, s *) 



X 


s 


where the condition ||a3;|| < 1, for all i £ Af, is dropped 
since the linearization is performed around x* = 0. Notice 
that 

14 ((B ® I)A(0, b*)(B ® I) T - D 2 ( 0, a*)) 1 3JV = 

= -1 t 3N D 2 (0,s*)1 3N 

= -353E(i- fl r®;)/: ( .»(2) 

l £ A/ - fc £ J\f i 

is zero if and only if s * = 1 N and negative otherwise (this 
includes all scenarios where at least two neighboring unit 


vectors are diametrically opposed), as long as /'(2) > 0 
for all k £ Ad. As such, the matrix in ( |XIII| i, which is 
symmetric, is not positive semi-definite for s * f 1 v (which 
means it has at least one eigenvalue with negative real part - 
Theorem 4.1.8 in [28]), which means any equilibrium where 
two neighboring unit vectors are diametrically opposed is 
unstable. On the other hand, the matrix in ( |XIII| > is positive 
semi-definite for s* = 1 N (this corresponds to the scenario 
where all unit vectors are synchronized), and the linearization 
cannot be used to infer stability of this equilibrium (nonethe¬ 
less, this equilibrium is asymptotically stable). 


XIV. Convergence to Constant Unit Vector 


Theorem [16] provides conditions for convergence to a 
synchronized network, but it does not provide any insight 
on whether lim n^f) exists, i.e., whether the synchronized 

t—f OO 

network converges to a constant unit vector or if it converges 
to a time-varying unit vector. 

Example 7: Consider the following unit vector, 
n (f) = [cos(ln(f))) sin(ln(f))) 0] T 

whose limit as time grows to infinity does not exist. Nonethe¬ 
less, its time derivative is given by 


sin(ln(f))) 


o 

1 

O 


cos(ln(f))) 

cos(ln(f))) 

= 

i o o 


sin (In (f))) 

0 


0 0 0_ 


0 


= s(j[0 0 l] J j n(t) = 5(w(f))n(f) 

where <*>(•) converges to 0. This means that convergence of 
<*>(■) to 0 does not guarantee convergence of its unit vector 
to a constant unit vector. 


Proposition 30: Consider two functions on [a, b], h 1 (x) 
and h 2 (x), such that 0 < h 1 (x) < h 2 (x) for all x £ 
[a,b\. Then, if f a h 2 {x)dx is convergent, so is f a h 1 (x)dx 
(comparison test) [29]. 

Proposition 31: Consider a function h(x) on [a, b\. If 
fl’ \h(x)\dx is convergent, then so is j'‘ h{x)dx [29], 

Proposition 32: Consider a unit vector n with angular 
velocity uj. If u> converges exponentially fast to zero, then 
n converges to a constant unit vector. 

Proof: Consider an arbitrary constant unit vec¬ 

tor e £ S 2 (for example, the standard basis vec¬ 
tors in R 3 ). Then e T n(f) is a scalar function of time, 
and it satisfies e T n(t) = f* e T S (u>(r)) n(r)dT. Since 
|e T 5 (w(t)) n(r)| < ||uj(t)|| < Cexp(—A t) (for some 

that lim e T n(f) exists. Finally, since e is an arbitrary unit 

t—>o O 

vector, lim n(i) also exists. ■ 

t—f OO 

The goal in this section is to study whether a network that 
converges to a synchronized configuration actually converges 
to a constant unit vector (in an unknown fixed orientation 


positive C and A), it follows from Propositions B0] and 31 






































frame). An analysis for the dynamic case (control at the 
torque level) is not performed, but rather an analysis for 
the kinematic case (control at the angular velocity level). 
Consider then that uj is the control variable chosen as uj 
K T (B® l)e. 

Consider that g k is of class V U V° for all k £ M. Also, 
assume that, for each edge k, there exists an interval I = 
[0 6], for some 0 < b < n, such that if dg( k n, E n) £ / then 

QL k ||e fc || 2 < d k { k n, - k n) < a k ||e fc || 2 , (XIV.l) 
where e k = g k (dg( k n, s n))«S ( fc n) - k n. 

Then, if the network converges to a synchronized configu¬ 
ration, i.e. rpl = ■ ■ ■ = n,-iV, it follows that lim n,(f) = n* 

£—>■00 

for some constant n* and for all i G N. 

In order to come to this conclusion, consider the Lyapunov 
function V(n) = d k ( k n, E n ) whose time derivative is 

rendered V = — ||(B ®I)e|| < 0 by the chosen control law. 
If the network converges to the equilibrium configuration 
where n, 1 = = n,A, it follows that V(t) converges 

asymptotically to zero, and consequently, for all edges k £ 
AL dg( k n, k n) will eventually enter the set I (which has 
non-zero measure). Say T is the time instant when, for all 
edges k £ M, dg{ k n, E n) £ I for all t > T. From the 
previous discussion, a finite T exists, and for t >T one can 
say 

M M 

^2a k ||e Jfc ( fc n, E n)|| 2 <V(n) < ^ a k ||e fc ( fc n, £ n)|| 2 

fc=l fe=l 

min(aj ||e(n)|| 2 <V(n) < max(a t ) ||e(n)|| 2 . 

k k 

Since V = — \\(B ® I)e|| 2 < —A m j n (B T B) ||e|| 2 , it follows 
that ||e|| converges exponentially fast to 0 (tree graph), and 
so does uj, (since uj = 7 Z T {B (g> I)e) for all i £ Af. As 
such, it follows from Proposition [32] that all n, converge to 
a constant unit vector. 


We can perform a very similar analysis for the dynamic 
case, but we need to construct a different Lyapunov function 
(than that provided in Section IV-D| >, specifically one whose 
time derivative depends on both || {B <g> I)e|| and ||. 
We do not present such Lyapunov function in this paper. 
Convergence to a constant vector in arbitrary graphs is also 
a topic for future research. 


We present next an example of a distance function that 
satisfies conditions (|XIV|). 


Example 8: Consider the distance function presented in 
Example [5] with a = 1, i.e., (((n^n,) = a(l — n)'n 2 ). 
For this distance function g = a and ||e(n 1 ,n 2 )|| = 

a ||S (rii) n 2 ||. This distance function satisfies 

1 ||e(n 1 ,n 2 )|| 2 , , 1 „ , ^ l|2 

-x- < d{n u n 2 ) < - ||e(n 1 ,n 2 )|| 

a Z a 

for n 2 ) £ [0, f]. 


XV. Extra Remarks 

Remark 33: In Definition [3] since /(•) is monotonous 
in (0,2), lim s ^ 0 + /(s) and lim s _ >2 - f(s) exist; while 
lim ( ,_ >0 + f'(s) and lim s _,. 2 - f'(s) do not necessarily ex¬ 
ist. However, in what follows the important limits are 
lim ( ,_ >0 + f(s)^/s and lim^,,- f'(s)y/2 — s which exist: 
lim 5 _, 0 + f'(s)yfs = 0 since lim sup s _> 0 +f'(s) < 00 and 
liminf^^o-i- f'(s) > 0; lim s _ >2 - f(s)\/ 2 — s = 0 when 
limsup s _ >2 - f'(s) < 00; and lim„_ ( , 2 - f'(s)\/2 — s may 
or may not be 0 when limsup s _ >2 - f'(s) = 00. Also, 
if / £ V°° then lim s _ >2 - f'(s)y /2 — s ^ 0: indeed, 

suppose lim s ^ 2 ^ f'(s)y /2 — s = 0, which implies that 
sup ag(0 2) /'(s)V2 — s =: b < 00; since 

lim f(s) = lim f(s) + lim f f'(x)dx 

•-*■2- —K) + Jo 

= lim [ J- - (/'( x)yj2 - x)dx 

=-*• 2 - Jo V 2 — x 

f s 1 

<b lim / d.x = 2v/2 b < 00 , 

~ s^-Jo sj2-x 

which implies that / ^ V°°. As such, / £ V°° => f £ 
V°. Figure [ 5 ] illustrates how the properties that /(•) satisfies 
affects the classes it belongs to. 

Definition 5: Given er : R 3 1 —> R 3 , we say er £ X if 
3(7 £ C 1 (R >0 ,R >0 ) : rr(x) = cr(||x||)x, and cr ,max = 
sup xeK 3 ||L>cr(x)|| = sup xeB3 ||^|| < 00 (i.e., cr(-) is 
Lipschitz, see Lemma 3.1 in [30]) Moreover, we denote 
^max _ SU PxeR 3 [|cr(x)[|. 

Remark 34: Theorems [16] and [24] provide asymptotic re¬ 
sults, namely linp^^ e(n(f)) = 0 (or alternatively that 
linit-joo dist(n(t), fi*) = 0, where fl* = {n £ (iS 2 )" : n = 
ljv (g> n*,Vn* £ <S 2 }). An open question is when can it be 
guaranteed that those limits converge to 0 exponentially fast. 
Similarly to [11], consider 9 X = f{\ — cos(0 1 — 0 2 )) sin(0 2 — 
9J) and 0 2 = /'(1 — cos(0! — 9 2 ))sin(9 2 — 9J), where we 
omitted the time dependencies; it follows that, if we denote 
9 = 9 1 — 9 2 , then 9 = — 2/'(l — cos(0)) sin(0). In [11], 
only /(s) = a 12 s => f'(s) = a 12 is considered, but in this 
manuscript we consider a wider set of functions. Consider 
then the function /(s) = |(y/s(2 — s)(s — 1) + arccos(l — 
s)), and notice that f £ V 0 (see Fig. [6j». It follows that 
9 = — | sin(0)| sin(0); therefore, if 0(0) £ [0, ^], it follows 
that 9{t) = arccot(f + cot(0(O))), which does not converge 
exponentially fast to 0 (since J o °° 9(t)dt = 00 ). Consider 
now the function f(s) = s, and notice that / £ 'Pg. It 
follows that 9 = —2 sin (9) < —9 , where the inequality 
follows for 9 £ [0, f ]; therefore, if 0(0) £ [0, |], it follows 
that 0(f) < exp(—f), which does converge exponentially 
fast to 0. This suggests that, in this manuscript’s framework, 
exponential convergence cannot be guaranteed for all / £ V, 
while asymptotic stability can still be guaranteed. 





















